THE QUARTERLY JOURNAL OF 


MECHANICS AND 
APPLIED 
MATHEMATICS 





Editoria’ Board 


Ss. GOLDSTEIN R. V. SOUTHWELL 
G. I. TAYLOR G. TEMPLE 


together with 


A. C. AITKEN; S. CUAPMAN; A. R. COLLAR; T. G. COWLING; 
Cc. G. DARWIN; W. J. DUNCAN; A. A. HALL; D. R. HARTREE; 
WILLIS JACKSON; H. JEFFREYS; J. E. LENNARD-JONES; M. J. 
LIGHTHILL; N. F. MOTT; W. G. PENNEY; A. G. 
PUGSLEY; L. ROSENHEAD; ALEXANDER THOM; 


A. H. WILSON; J. R. WOMERSLEY 


Executive Editors 
G. C. McVITTIE Vv. C. A. FERRARO 





VOLUME IV 
1951 


OXFORD 
* THE CLARENDON PRESS 





Oxford University Press, Amen House, London E.C. 4 


GLASGOW NEW YORK TORONTO MELBOURNE WELLINGTON 
BOMBAY CALCUTTA MADRAS CAPE TOWN 


Geoffrey Cumberlege, Publisher to the University 


PRINTED IN GREAT BRITAIN 
AT THE UNIVERSITY PRESS, OXFORD 
BY CHARLES BATEY, PRINTER TO THE UNIVERSITY 














THE QUARTERLY JOURNAL OF 


MECHANICS AND 
APPLIED 
MATHEMATICS 


VOLUME IV PART 1 
MARCH 1951 


UNIVERSITY 
OF MICHIGAN 


MAY 4.1951 


ENGINEERING 
LIBRARY 


OXFORD 


AT THE CLARENDON PRESS 
1951 


Price 12s. 6d. net 


T BRITAIN BY CHARLES BATEY AT THE UNIVERSITY PRESS, OXFORD 





THE QUARTERLY JOURNAL OF MECHANICS 
AND APPLIED MATHEMATICS 


Editorial Board 


S. GOLDSTEIN R. V. SOUTHWELL 
G.I. TAYLOR G. TEMPLE 

together with 
A. C. AITKEN J. E. LENNARD-JONES 
S. CHAPMAN M. J.. LIGHTHILL 
A. R. COLLAR N. F. MOTT 
T. G. COWLING W. G. PENNEY 
Cc. G. DARWIN A. G. PUGSLEY 
W. J. DUNCAN L. ROSENHEAD 
A. A. HALL ALEXANDER THOM 
D. R. HARTREE A. H. WILSON 
WILLIS JACKSON J. R. WOMERSLEY 
H. JEFFREYS 


Executive Editors 
G. C. McVITTIE Vv. C. A. FERRARO 


THE QUARTERLY JOURNAL OF MECHANICS AND APPLIED MATHEMATICS is 
published at 12s. 6d. net for a single number with an annual subscription 
(for four numbers) of 40s. post free. 


NOTICE TO CONTRIBUTORS 


1. Communication, Papers should be communicated to one or other of the Executive 
Editors, by name, at King’s College, Strand, London, W.C. 2. 


2. Presentation. Manuscripts should preferably be typewritten, and each paper should 
be preceded by a summary not exceeding 300 words in length. References to literature 
should be given in standard order, author, title of journal, volume number, date, page. 
These should be placed at the end of the paper and arranged according to the ondedd of 
reference in the paper. 


3. Diagrams. The number of diagrams should be kept to the minimum consistent with 
clarity. The lines of the figures should be drawn in ink either on draughtsman’s paper 
or on good quality white paper. Each individual line in the figure should bear reducing 
to one-half of the size of the original, and great care should be exercised to see that the 
lines are regular in thickness, especially where they meet. Lettering of the figure should 
be in pencil and should be sufficient to define clearly the lines and curves in it. The 
writing of formulae or of explanations on the diagram itself should be avoided. All 
explanations of symbols, etc., should be given in underline. Contributors should indi- 
cate on their manuscripts where figures should be inserted. 


4. Tables. Tables should preferably be arranged so that they can be printed with the 
columns parallel to the longer edge of the page. 


5. Notation. All single letters used to denote vectors in the manuscript should 
be marked by underlining with a wavy line. Scalar and vector products should be 
denoted by a.) and a » b respectively. Real and imaginary parts of complex quantities 
should be denoted by re and im respectively. 


6 Offprints. Authors of papers will be entitled to 25 free offprints. This number is 
available for sharing between authors of joint papers. 


7. All ees eorygumed other than that dealing with contributions should be addressed 
to the Publisher: 
GEOFFREY CUMBERLEGE 
OXFORD UNIVERSITY PRESS 
AMEN HOUSE, LONDON, E.C, 4 








TI 
bour 
on t 
it is 
the 
nam 
first 


com 


for 


touc 
tion 
mo! 
and 
0-09 
valt 
thar 
den 
the 
T 


fron 


hot p 





THE LIFT AND MOMENT ACTING ON A 
CIRCULAR-ARC AEROFOIL IN A STREAM 
BOUNDED BY A PLANE WALL* 


by 8S. TOMOTIKA, K. TAMADA, and H. UMEMOTO 
(Physical Institute, Faculty of Science, University of Kyoto, Kyoto, Japan) 


Received 25 July 1950| 


SUMMARY 
The lift and moment of an aerofoil in the form of a circular are placed in a stream 
bounded by a plane wall are discussed. The earlier work of Green (1) and Hudimoto (7) 
on this problem is reviewed and by adopting a method similar to that of Hudimoto 
it is found that when the camber of the aerofoil is small the effect of the wall on 
the lift and moment is similar to that for an aerofoil in the form of a flat plate, 
namely, as the arc-aerofoil approaches the wall the lift and moment coefficients 
first decrease and then increase to values which are greater than the corresponding 

values for an are-aerofoil in an unlimited stream. 


1. Introduction 
SoME years ago. Green (1) investigated the manner in which the ground 
effect upon the lift and moment of an aerofoil is modified by its camber, 
by evaluating the forces acting on a circular-are aerofoil placed in a semi- 
infinitely extended stream in any position near the bounding plane wall. 
He performed a few numerical calculations for a particular circular-are 
aerofoil with camber approximately equal to 0-097 in the comparatively 
simple cases when the centre of the circular-are aerofoil lies on the bound- 
ing wall and when the circle on which the are lies touches the wall. By 
combining his results with the previous results obtained by one of us (2, 3) 
for a limiting case in which the trailing edge of the circular-are aerofoil 
touches the wall. Green has indicated the general nature of the modifica- 
tion introduced by the camber on the ground effect upon the lift and 
moment forces acting on an aerofoil. He was led to suppose that the lift 
and moment coefficients of a circular-are aerofoil (with camber equal to 
0097) are probably always decreased by the presence of the wall for all 
values of the ratio / H and for all values of the angle of incidence « less 
than 85° (including the small values which oceur in practice), where / 
denotes the chord of the aerofoil and H the distance of the midpoint of 
the chord from the wall. 

These results differ from those for a flat plate (plane aerofoil) and also 
from experiment (except when the distance of the aerofoil from the ground 


+ This paper was « pleted by the end of March 1942. Publication at that time was 
hot possible, howe 
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is large). In previous papers (4, 5, 6), one of the present writers has shown 
that for small values of the angle of incidence « (e.g. « = 5°, 10°), as the 
flat plate approaches the wall, the lift coefficient first decreases and then 
increases to values which are greater than the values for a flat plate in an 
unlimited stream and that similar results also hold for the moment 
coefficient. 

A flat plate is evidently the limiting form of a circular-are aerofoil 
whose camber tends to zero. It seems therefore probable that the results 
for an aerofoil in the form of a circular are with a camber smaller than that 
considered by Green will differ from his results, but will be somewhat 
similar to those for a flat plate. Thus, in order to determine the true nature 
of the effect of the wall upon the lift and moment acting on a circular-are 
aerofoil, it is desirable to extend the investigation to cases of aerofoils with 
fairly small camber; and the principal object of the present paper is to 
present the results of our detailed numerical calculations. 

We have reinvestigated the problem ab initio by employing suitable 
conformal transformations, which are somewhat different from those used 
by Green. The analysis we have developed is similar to that given ina paper 
by Hudimoto (7), who has also discussed the effect of the ground upon 
the lift and moment of a circular-are aerofoil. His results for one special 
casey are similar to Green’s, namely that the lift coefficient of the circular- 
are aerofoil is decreased by the presence of the ground. 

The final expressions obtained by Hudimoto for the lift and moment are 
unfortunately inelegant and complicated. However, we have been able to 
obtain the general formulae for the lift and moment in fairly elegant forms 
as in Green’s paper, although we have employed a method of analysis 
similar to Hudimoto’s. 

We have carried out detailed numerical calculations for three circular- 
are aerofoils, the cambers of which are respectively 0-022, 0-053, and 0-097 
approximately, the value of the angle of incidence « being taken to be 5 
in all cases, and for values of the ratio //H between 0 and 22-95. 

Thus it has been found that when the camber is sufficiently small and 
the angle of incidence is also fairly small, as in practice, the manner in 
which the lift and moment coefficients of a circular-are aerofoil are affected 
by the presence of the wall is similar to that for the case of a flat plate: 

namely, as the circular-are aerofoil approaches the wall, the lift coefficient 
first decreases and then increases to values which are greater than the 
In his numerical calculations a circular-are aerofoil with camber equal to 0-115 was 
placed at an angle of incidence 17° 2-28’ and the value of the ratio //H was 0-885. 
No numerical discussions on the moment have been made by Hudimoto. 


§ The latter value 22-95 corresponds to the limiting case when the trailing edge of the 
are touches the wall. 








valu 


hold 


2 3 

fie 
as tl 
invis 


an i 


aero 
(Fig 

Tl 
dete 
Hud 
of tl 
HH’ 

D 
tO al 


The 
of bi 


coin 


Furt 

den 
N 

we: 

















LIFT AND MOMENT ACTING ON A CIRCULAR-ARC AEROFOIL 3 


values for the circular-are aerofoil in an unlimited stream. Similar results 


hold also for the moment coefficient. 


2. The conformal transformations 

Taking the plane of the fluid motion (assumed to be two-dimensional) 
is the z-plane, we consider a steady irrotational flow of an incompressible 
inviscid fluid past a circular-are aerofoil, AA’, placed in any position near 
in infinite plane wall HH’, and such that the circle on which the are 





ly 
A 
J / 7 
' ' 
\ a a I Bi 
H* X O . 
Fic. | 


verofoil AA’ lies intersects the bounding wall at two real points X and Y 
Fig. 1 

The complex velocity potential for the flow under discussion will be 
determined by using certain conformal transformations introduced by 
Hudimoto. We take the origin O of the z-plane to be the midpoint of the join 
of the intersections X and ), and we take the real 2-axis along the wall 


HH’, while the arc-aerofoil is assumed to lie in the upper half of the z-plane. 





Denoting by a the length of OX and OY, we transform the z-plane on 
to an f plane bv the transformation 
2 a 


f log ; (1) 


Z a 


(he upper half of the z-plane is then transformed into an infinitely long strip 
if breadth 7 parallel to the real axis of the f-plane, one side of the strip 
oinciding with the real axis of the f-plane. The points X and Y in the 
:-plane correspond to the point at infinity in the f-plane, while the point 

infinity in the z-plane is transformed into the origin of the f-plane. 
Further. the circular are AA’ is transformed into a straight line-segment 
denoted as AA’ parallel to the real axis. The f-plane is shown in Fig. 2. 

Next by making a cut along CGG’'C’. as shown in Fig. 2. where both 


CG and C’G’ pass through the midpoint ./ of AA’ and are perpendicular 
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to the real axis, we transform the f-plane on to the upper half of a t-plane 
(Fig. 3) by the method of Schwarz—Christoffel. The transformation is 
if ae 
a M 9 9 9 ‘ 9 9 9 ? (2) 
dt (t?—h?),/{(t?—c?)(t2 —g?)\’ 











: cele HH (f=0) ae 4 


lic. 2. f-plane. 


| 
| 
: 
| 
RY GCA M ACG XR 


ric. 3. ¢-plane. 
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Fic. 4. s-plane. 


where WV is a certain constant which will be determined presently, while 

b,c. g. h are all reai constants and the points t = b, —b. ec, —c. g. —q. h. 

and —h correspond to A’, A, C’. C, G’, G, X, and Y respectively. 
Lastly. by introducing Weierstrass’s (9 function with periods 2w,. 2ws, 


where w, and w,7 are real and positive. we conformally transform the 
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upper half of the t-plane into a rectangle of sides 2m, and ws/7 in an s-plane 


by the relation 


{2 (2(s)—es (3) 
The points A, A’, C, C’, G. G, X,Y, M, R now correspond to s m, — Bp, 
wW,+Ws W,+W3, WW, —W, v, Vv, Wg, 0 respectively, while both the 
points H. H’ correspond to the same point s = 89. say (Fig. 4). 


3. Now. from (2) and (3) we have 


af M (Ss) (2) 
ds —(s)—O(r 


y PW) () (2) 


\ 


[C(u+v)—C(p—v) —C(s +r) + C(s—v)]. (4) 
{J (Vv) 


This differential equation can be integrated to give 


f us | C( py v) C(p v)is log & by . (5) 


Nv) (pL) | ce a(s+y) 

()'(v) a(s v) 
where C' is an integration constant to be determined presently, and o(s) 
is the Weierstrass sigma function. 

Since the function f is continuous along the cut CGG'C’ artificially made 
in Fig. 2. it must necessarily have a period 2w,:; thus the function f must 
satisfy the relation f(s+-2w,) = f(s), and this condition gives immediately 
. relation between the constants p and v in the form: 


u-+-1 C(u—v) 291 V/@}. (6) 
where, as usual, », Wy 
\oain. since the values of f at R (s 0) and at G (s w,) differ by iz. 
we have i 
f, f(w,)—f(0) (df ds) ds iT. 
0 


Thus. integrating the function df/ds from s Oto s w,, bearing in mind 
that the point s vis a simple pole of the function df ds and taking 


6) into account. we get 


Q'(v) 


Vl R 
V(v)— Pp) 


(7) 
Combining this result with (6), the expression for f can be written in the 
lorm ( ) 
ot O\s V 
1” s—log C. (8) 
, a(s v) 
The constant C on the right-hand side can be determined by the con 
lition that the points H, H’ correspond to f 0 in the f-plane and to 
‘) in the s-plane. We thus have 
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and on substituting this value of C' in (8). we have ultimately the expression 
for f in the form 


> 2,1 o(s+v)a(s v) 
J 2 So) log ° 
Wy o(s—v)o(S)+v) 
log U1els v) wf {4(s6 + py) W1} (9) 
Byia(s+yr) 3711 3(Sp v) wi} 


4. Denote by @ the angle which the chord XY subtends at any point 
of the are AA’ in Fig. 1. Then the complex coordinate of the point A’ in 
the f-plane may be written as f, p+70, where p denotes the distance of 
A’ from the imaginary axis of the f-plane (see Fig. 2). Also, we denote by 
2d the distance between the two points A and A’. These three quantities 
6, p, and d can be expressed in terms of , v, and s, as follows. 

Firstly, since the values of f at C (s w,+ws) and at G (s w,) differ 
by i@, we have 

fo—fe f(w,+a3)—f(w,) 10. 
and when use is made of the second expression for f in (9), this condition 


sives immediately 
¥ . : 0 TV / Wy. (10) 


Next. since the distance between the points A and A’ is equal to 2d and 
M is the midpoint of the sé¢gment AA’, the values of f at s —= ws and at 
s = p evidently differ by d. We therefore have 

fu-—Si f(ws)—f (uv) d. 
and this condition gives 


3, (37 —bv/w,)I fa (u+v)/o} (11) 
d (47+ bv/w,)9 4 3(u—v) w,) 
where, as usual, T = w3/a. 

This expression for d can be greatly simplified, however, if we make use 
of the well-known properties of theta functions. In fact, putting 


ph m—- We. ( ] 2) 
where m is real. we have 
951 ) 
PAt(m+v)/w 
d log ai al 1s 


© ¢ fl es (13) 
Fi d(m v)/@4; 


Also, bearing in mind that the value of f at the point A’. which corre 
sponds to s ji in the s-plane, is given by f, p+ié we have from 
(10), (12). and (13). 


} +V)/ Wy} 


p+ié log 


Bf dr—F(m+v)/w,}F 44 (s9 
h{dr—3(m v)/wyjP1{3(S9—v)/@,} 


Cc 
FS (spt+v)/w,} 


v)/wy} 
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from which we obtain 


\ fl > ' 
V4, 9(SgV)/ 45 


p—d log (14) 


Hy) 3(Sy—v)/04} 
Further. for convenience we shall rewrite the relation (6) in terms of 


theta functions. using (12). We then have 


# f1(m+v)/w,! # {i(m—v)/w,! of 
12 1) 412 1 (15) 


aA3(n V)/ Wy) I d(m v)/@45 


, 


From the given geometrical properties of the circular-are aerofoil AA 
we first calculate the values of the three quantities 6, p, and d, using 


certain relations which connect these three quantities with geometrical 





parameters defining the circular arc. This will be discussed subsequently. 
lhe values of the three quantities m, v, and s, can then be found by solving 
equations (10), (13), and (14), while equation (15) determines the value of 
the parameter g (= exp(rz?)) for the elliptic functions as well as for the 


theta functions used here 
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5. Geometrical parameters defining the circular-arc aerofoil 

We are now in a position to consider the geometrical parameters defining 
the circular-are aerofoil AA’. Referring to Fig. 5. let | be the length of the 
chord AA‘, R the radius of the are AA’, 2y the angle which the are AA’ 


| subtends : 


t its centre O’, and H the distance of the midpoint of the chord 
1A’ from the bounding wall. Further, let « be the acute angle which the 
chord AA’ of the aerofoil makes with the wall. Then « is the angle of 
incidence of the aerofoil when the fluid flows in the positive direction 

| of the x-axis. Since 6 denotes the angle which the chord XY subtends at 

iny point on the are AA’, both the angles XO’O and YO’O are evidently 

equal to 6 


Further. remembering that OX OY a. we have the following 


oby ious relat ions 


Rsin@ =a Rsiny = W. 
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from which we get 
sin @ a l 
a pean teaeh R ; ‘aT 
“ siny sin 6 “ siny 
Also, it is easily seen that 


H R(cos 6+-cos «cos y), 


and the maximum height of the circular are AA’ is equal to R(1 


UMEMOTO 


(16) 


(17) 


Cos y). 


Thus. the so-called camber o of the arc-aerofoil, which is defined as 


usual by the ratio of its maximum height to the chord length /, is given by 


R(1—cosy) 1—cosy 


o ——* i tan dy, 
l 2siny Z . 
and the ratio //H becomes 
l 2siny 
H cos 6+ cos «cosy 


/ 


(18) 


(19) 


Next, consider the expression for the distance 2d between the two points 


A and A’ in the f-plane. Let the complex coordinates of these points be 


denoted, as before, by f, and f, respectively, and let those of the corre 


sponding points in the z-plane be denoted by z, and z, respectively. 


Then, bearing in mind that the f- and z-planes are connected by the 


relation (1), we have 
. . Z A)(z a 
2d bav—Ia log 4 Nz ) 
’ (z,—a)(z4+a) 
7 4 
and after some reduction we have ultimately 


COS a cos(@ y) 


d log . 
~ cos «+ cos(@+-y) 
Further, since f, p+ié, we have 
. ‘ yA a P 
p fy— log 26, 
Z a 
1 
z,—a)(z a : 
and therefore p—d 1 log —4 Nzq 10. 


(Z4+@)(Z4-+@) 


After some reductions we get the expression for p—d in the form: 


ee walt 
p—d fcr A cos(@+- «) 


cos y+-cos(6— «) 


Combining (20) and (21) with (13) and (14) we have ultimately 


Hi d(m tp) W1} {COS a> cos(@ y)\3 
44 3(m—v)/a} \cos «+cos(4+-y)} 
{Msp t+v)/ar4} _ (c0sy-+eos(0-+a))J 


Fi{3 (Sp v) Wy} \cos y ! eos(A—«)} 


(20) 


(22) 
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and if the values of «, y, and //H (or @) are given, then these two equations, 


together with equations (10) and (15), determine the values of m. v, s5, and q. 


6. The complex velocity potential 

We shall now obtain the complex velocity potential for a steady irrota- 
tional continuous flow of an incompressible inviscid fluid past the circular 
are aerofoil AA’ placed in the neighbourhood of the wall HH’ (Fig. 1). We 


issume that the fluid at infinity flows with constant velocity U’ from left 











to right. i.e. in the positive direction of the 2-axis. Also, we assume for the 
present that the circulation round the aerofoil is zero. Then the flow 
pattern will become as shown schematically in Fig. 6. 

Let w, be the complex velocity potential for the flow under considera 
tion. Then, in the z-plane the conjugate complex velocity dw, /dz is every- 
where finite except at the two edges A and A’ of the arc, where, as will 


be seen later. it becomes infinite. 


If now we transform the z-plane on to the f-plane (Fig. 2), then as we 
have seen the point 2 © is transformed to the origin of the f-plane. and 
mn taking account of the obvious relation 

du, df (dw, dz)(dz df ). 
together with an expansion 
dz/df 2a f*?+ positive integral powers of f. 


which follows immediately from (1) on assuming f to be very small) it is 
readily seen that the origin f — 0 becomes a singular point of the flow in 
the f-plane 

On the other hand, the two poles at the leading and trailing edges of the 
we-aerofoil are transformed respectively to the ends A’ and A of the line- 
segment AA’ of the f-plane, and it can easily be shown that these points 
are still the poles of the function dw,/df. 
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Also, it is readily seen from (1) that (dz/df), 
also that f = 0 corresponds to the points X (z —a) and Y (z = a) in 


cq = 9, and remembering 
the z-plane. we have (dw,/df), , 0. The flow pattern in the f-plane is 
shown schematically in Fig. 7. 

If we further transform the f-plane on to the s-plane, the pole at 
f = 0 (H.H’) is transformed to the point s = sy. However, since it is 
readily seen from (4) that df/ds = 0 when s = pand s ji. the function 


dw,/ds has no longer poles at the points A and A’, but is rather finite and 
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regular there. Furthermore. in spite of the fact that (dw, df), , 0. 
dw,/ds is finite and regular at the points X and Y. since (df ds), ,, x 


as will be seen from (4). 

Thus the conjugate complex velocity dw, /ds in the s-plane is everywhere 
finite inside the rectangle except at the point s = s, where it has a pole 
of certain order. 

We shall now consider the order of a single pole at s = s, of the function 
dw,/ds. First. we have 

dw, dw,dzdf dw, dz(df\? ds 

ds dz df ds dz df\ds} df 
Since the velocity of the undisturbed uniform flow at infinity in the 
z-plane has been denoted by l’, we have (dw,/dz),_, U. Also, as shown 
previously, we have 

dz/df — 2a f*+-positive integral powers of f, 
while equation (9) gives immediately 
f x(s—s,)+B(s—sy)?+... (x #0). 


Therefore. remembering that z = 0 corresponds to f = 0 as well as to 


dz(df\* 2a(' df\? 2a 
i a ae 


3 = 8. we have 
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und consequently, when s tends to s 


0° 
dw, 2aU es a £ ( ) (23) 
OSItIVe Integra owers O Ss Sy). yey) 
ds (s—s,)? dj). I ] 0 
Since, however, (ds df), . is finite, we find ultimately that the function 
lw, ds has a pole ot the second order at s ‘o- 


Further, dw, ds must be real on both the lower and upper sides (@G’ and 
CC’) of the rectangle in the s-plane, which correspond respectively to the 
rigid plane wall and the are aerofoil in the z-plane. Hence, it is seen, by 
Schwarz’s principle of reflection, that dw,/ds must have a period 2ws. 
Moreover, it is readily seen that the said function dw,/ds must have one 
nore period 2w,, since it takes the same values at any corresponding 
points on the two sides CG and C’G’. 


Thus we find that dw, ds is a doubly-periodie function with periods 2w,, 


2w, having a pole of the second order at the point s = sy. Such a function 
in, however, be uniquely expressed by using Weierstrass’s ~ function 
with periods 2w,, 2w, having a pole at s = s. We thus have 

dw, ds uw. (Is So) + C1; (24) 


vhere w, and (, are constants which can be determined as follows. 
Firstly. the constant (, can be determined from the condition of no 


irculation round the aerofoil. which can be expressed in the form 


(dw, ds) ds (). 


) 


Inserting the expression (24) for dw,/ds in the integrand and performing 


the integration we have 
Cy Ws 11/4. 


lhus the expression for dw, ds becomes 


dw, ds W.4(S— Sq) + 74/4}. (25) 

Next. the constant w, can be determined by comparing the expansion 
23) for dw, ds with that for (25). From (25) we have 

ds s—s,)* 


0/ 


positive integral powers of (s— sp). (26) 


nd hence it is easily found that 


u Yal'(ds/df). .. (27) 
However. we have from (8) 
lf = 2,1 , 
: vj)+ O(s—v) 
is Wy 
| V3i{9(S—v)/wy, t 4 1(s-yp) Wj 


v) Wy} 
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Thus the expression for w, becomes ultimately 


9’ 4 


as 
Fid(Sp—v)/@,$ Ayd(Spt+v)/a,}] 


Ww. dal wy - a5 ; (28) 
Fyd(Sp—v)/a,}  Ayd(Spt+v)/a,} 


The flow pattern in the s-plane is shown schematically in Fig. 8. 
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7. From the relation 

dw dw, ds df 

dz ds df dz 
and by taking into account the results that ds/df is infinite at s — p and at 
s ji = 2wg—p. while dw,/ds is in general finite there, and that df/d: 
is also finite at the corresponding points A and A’, we find that for the 
flow with complex potential w, so far discussed, and for which there is no 
circulation round the arc-aerofoil AA’, the fluid velocity becomes infinite 
at its two ends A and A’. In order to avoid the infinite fluid velocity at 
the trailing edge A of the aerofoil, we superpose on the above flow a 
rotational flow of circulation « round the circular are, choosing « so that. 
as usual, the fluid velocity at the trailing edge becomes finite. The fluid 
will then leave the trailing edge smoothly. 

We assume that the circulation « round the are-aerofoil in the z-plane 
takes place in the clockwise sense, as would be the case if the aerofoil is 
placed in the stream as shown in Fig. |. It is evident that this circulating 
flow becomes, when transformed to the s-plane, a uniform flow parallel to 
the negative direction of the real axis. The complex velocity potential for 
this circulating flow. which satisfies the boundary conditions, is therefore 


given by 


We os Ss, (29) 
The corresponding conjugate complex velocity becomes 
dw, ds Cas (30) 


where C, is a positive constant. 


In the general case when there is circulation round the aerofoil. the 
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complex velocity potential for the continuous flow past the aerofoil is 
obtained by adding the complex velocity potential w, to w,. Denoting 


this by w. we have , , 
this by w W,+ Ws w,—C,8, (31) 


und so. from (25) and (30 we get 


dw/ds W4G(S— 8) + Hy /@443—Co. (32) 

The constant C, can be determined by the Kutta condition that the 
fluid leaves the trailing edge of the arc-aerofoil smoothly. Bearing in mind 
also that the trailing edge A corresponds to s = pw in the s-plane, the 


Kutta condition may be written as (dw/ds),_,, = 0, and therefore we have 


(> Wi (Le — 89) + 9y/@45- (33) 


If we substitute this in (32), we have ultimately 
dw/ds W§9(S— Sy) — (t— Sp) j- (34) 


In this case the circulation « becomes, by (32), 
K | (dw ds) ds 2w, ( 2 2w, Wy (pe So) +1, Wy}; 


ind if use is made of theta functions. we have. after some reduction. 


u oS lin Sp )/ Wy; ie 
« |Mate wd id, (35) 
Ft d(m So) Ww} 





HAd(m So)/@4} *) 

2w, lA d(m Sq) / Wy; | I 
where m is the real quantity previously defined in (12) such that p = m+-ws. 
Further, if we use the value of w, given by (28) together with the first 


relation in (16). we have 


te 


al Sin 4 | hy d(m S,) Ww Fi 3(m So) W4} 


0 


SIN 6 (4h 3(M— so) wy, | Paid (m small 
9/51 , ' 9’ $1 i i 1 
{Vy 9(Sy—V)/ 4; V9 (S97 v)/wy5| 


« 1 A 
LA13 (Sp V)/ wy) FyLa (Sp rv) w}) 


8. The general expressions for the lift and moment forces acting 
on the circular-arc aerofoil 
If. as is usual for a very thin aerofoil, we do not take into account the 
singularity at the trailing edge. the lift on the cireular-are aerofoil can be 


evaluated as usual by Blasius’s first formula. Denoting the 2- and y 


mponents of the resultant force acting on the aerofoil by X and Y 
espectively We have 
. [ (dw? ” 
X ) tp p | dz. (37) 
* dz 
{ 
vhere the integral is taken in the counterclockwise sense along anv closed 


ontour C surrounding the aerofoil and p is the density of the fluid. 



































14 S. TOMOTIKA, K. TAMADA, AND H. UMEMOTO 


The evaluation of the above integral can be conveniently carried out by 
transforming the integrand so that the integration takes place in the 
s-plane. Since, however, our calculations are similar to those in Green’s 
paper (1), only the final results will be given here for brevity. Thus, putting 








f S > 
= fs 0 bs (38) 
=W 27 =W, <7 
ve have 
. y a ee ee ee 
x 0. = —— C3 E*(Y,+Y,+Y3;+Y,+-Y;), (39) 
tol 2) sin y i , 
where 
y (FP(O)IES(P1+42)/7} , 93(0)) 
| 92(0)93{4(4,+¢4,)/7} | 3,(0)) 
i H2(0)9354(h,+¢5)/7} , H(0) —-FY(0) 
| | 93(0)9354(6,+¢5)/7} © 9,(0) — 394(0) 
, go { ds ar! 1 | 2) a! 
HP(O)| - (P1—P2)/ 7} o gat $a) | (40) 
9411 (6, —d)/77} wan 3($,+¢5)/77} 
, » pF (OVI f(b, + b9)/7} (84 {(b,+-40)/77} 2 Pali ¢,-+ os )/a') (41) 
- ie IE 1(¢,+¢5)/7} LA, + hg)/ 7} F413(b1 +) my 
: 3" d 17) ne by 7 
Y, {Va 4 Vi\P2 ) (42) 
\it,(d, Tr ) Ii (de a)} 
: 27 sin( y 
y, D ‘—y) (43) 
cos 6+-cos(a—y) 
y 7 cot 6 | #7(0) 5 (, +8) / 7} VEE >( $y O)/77 
C | HO) | FH4(d,+-6)/7}  FH3(b,—O)/7} 
TE (4, $2) Ts 5 (dy I 2) 7 
FF} (dy dy) TT} 3 5a 4 ( (4, T ps) 7 
ie[ pu 7) ut >(P4 + va (44) 
. (dy, 77) UE (4, Ts | 
c+) RU, — 0) (45) 
F?(O)F, (A, 7)P, (0/7 
Q/2 C ae Pee sites 2 
. #2 (0)9,(4, om (h,/77) (46) 
Pr 4 (6, T po) 7 ja (dy dy) w 
and 
BE F(0)9(b1/7) 913 (b2+9)/ tA 3(h2—9)/ 7} (47) 


8,64 1(h,+4)/7} P43 4(¢,—9)/7 aid 4(h,+45)/7}4{3(b,—¢b0)/75 


9. By means of Blasius’s second formula the moment I, of the forces 
acting on the aerofoil about the origin of the coordinate-axes, reckoned 
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positive in the clockwise sense, is given by 




















\ 
, , * (dw\? 
- re bp g | z dz, (48) 
es salts dz 
( 
e where, as before. the integral is taken in the counterclockwise sense along 
: the closed contour C surrounding the aerofoil. As in the case of the lift. 
| the evaluation of the above integral can be conveniently carried out by 
transforming the integrand so that the integration takes place in the 
()) -plane. Omitting intermediate analysis, the final result is 
‘i — sin?4 
— 7C* E?D -—.— : (49) 
LoU?l? sin“y,cos @+-cos(a—y)} 
where (', D. and E are given by (45), (46), and (47) respectively. 
Further, it is readily found that the moment [ of the forces about the 
trailing edge of the circular-are aerofoil is given by 
Y= 14+ Yz,, (50) 
0) where Lg Rsin(a y) Sl sin(« } y) sin y. (51) 
\ Oe 
| 
| 
| 
12 a X y 
Green's s-plan: 
13 i 
| | 
| 
' 
ae O ae ens ed 5 
Y(s) 0) X(s3) 1: Y(sy+2e,) 
N(2. —) A(fH) 
pa ------- ----- 
| 
+4 | our s-plane 
5 | 
| (0) 
r Hie Xe” Yo 
Fic. 9. 
10. Comparison with Green’s results 
17 We shall now compare oul expressions for the lift and moment acting 
on the aerofoil with the corresponding expressions obtained by Green (1). 
rces The correspondence between various points in the s-plane in Green's 


ned analysis and in our s-plane is shown in Fig. 9 and in Table I. Thus it is 
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easily found that Green’s @,. 6. 45. and @, correspond to our ¢,+-¢,—7. 
d,+¢.—7, $,.—0+7, and ¢,+6—7 respectively, where ¢,/7 = m/w. 
b/ 7 8,/a, O/a V/W)- 
Further, we readily find that Green’s 6 is equal to our 7—@ and that 
Green’s « (= $(@,—4,)) and ¢ (= 0,—6 = 6,+-8) correspond to our ¢, and 
d, respectively. The correspondence between various angles is shown in 


Table II. 


mn fa bl 
PABLE | PABLE II 
Green's Present Green's Present 
analysis analysis analysis analysis 
8) p—8&)— 2a, 7—0o 6 
Sy 2ws p So x d, 
o5 vate b $s 
Sy 1 89 — 2a, 


Taking these results into account, and bearing in mind that in Green’s 
analysis 6, r+d—z7 and @, .+6—7, we then compare our various 
expressions with those of Green. After some reductions, we found that. 
except for the difference in notation, our expressions for the lift and 
moment acting on the aerofoil are in perfect agreement with the corre- 


sponding ones of Green.t 


11. Special cases 

(a) Centre of the circular arc on the bounding wall 

The general formulae for the lift and moment forces acting on the 
circular-are aerofoil are very complicated, but they assume simple forms 
when the centre of the circular are lies on the bounding wall and when 
the circle on which the arc lies touches the bounding wall. 

We begin with the first case. i.e. when the centre of the circular arc lies 
on the wall. Then @ = $7 and therefore we have 

Vv Sw). ph yw Wes. m Say. (52) 


The first equation in (22) reduces to 


t.(0) cos «+ sin y)4 = 
= - < \2 (9) 
P,(0) \eos x sin y| 
Also, remembering that we may write }s)/w,+ } lejn (3 > (), 
we have. from the second equation in (22), 
#, (1/77) cos y—sin x)2 . 
a 4 { / ; \» (54) 
,(4€ 7) (cos y+sin a| 
while the expression (19) for the ratio /}H reduces to 
iH 2 tan y/Cos «a. (55) 
[t may be remarked here that there is a misprint in Green’s expression for Y,: the last 


term should be read as: #,(0,/7)/34(40,/7) instead of 3.(0,/7)/34(40, 7). 
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Finally, after a little reduction, the expression (39) for the lift coefficient 


pecomes } 
Lol D,(Y; T Y,+ Y34 Y,), (56) 
re 8) 
whe re 
P(Le/a)P5(4e/a ~ 
D I a = : ” (57) 
79S (0)99(L€/7)85(fe/77)sin y 
72 (O)9, (Le/7)3,( 4/77) _ 
E * a, (58) 
U(4€/ 7) 4(4€/7) 
(#7 (0)99(4e/7) H5(0)) 
| P(O)9(Le/7) © 3,(0)f 
(4p [OPI a) HO) (0) 
| | H(O)R(e/z) © 80) 394(0) 
‘ ] l | 
Q'3 ¢ 1 ~¢ 
Uy (OO) (s€/ 77) | = + : (99) 
, — H(te 77 ) F4(4e am) | I 
BY j 1 = GQ’ /1 a GQ’ (1 = 
} 2 BVI ) vile 7) {0 (s€/7) 9 Valse ™)\ (60) 
B(ie/z) =| 8, (Ge/z 4(1Le/77)} 
H4(0)92(0)9,(4e/7) 95 (Lela , 
Ys = ~ §i(de/m)94(de/m) (61) 
U'a( 4 €/ 77 )U'4( 4 €/ 77) 
Iau 4 ())92(0)92 (4 TT ) 
\ <7 2 2\o€/7), | ol gS 
ly 4 94 tan(a—y). (62) 


1 4\ té 7) 
\lso. the expression (49) for the moment [, about the origin of the 
rdinate-axes reduces to 
| H | DX, (63) 
Lol 2 sin ysin(a—y) 
is readily found that. except for the difference in notation, these 
xpressions for the lift and moment are in agreement with the corre- 


sponding expressions given by Green. 


The ( Le on which the are les touches the bounding wall 
In the second case in which the circle on which the are lies touches the 
inding wall. the quantities a and @ tend to zero in such a manner that 
R. where R is the radius of the circular are. Since 6 = zv/w,. v tends 
zero also and equation (15) reduces to 
9 mw.))? H (1m w,) ' 
| a\° l | - 2 1 ; (64) 
[U4 yt Wy) |} H4(3m W,) 
is before. we put m Q,/7 and So Wy, hy 7, this equation 
COmMmeS 
: (65) 


« r/* ; F4(24, 7) 


 ( id, Tr ) 
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while the limiting forms of the two equations in (22) become respectively 


H,(4,/7) 7siny ~~ 
a71 = 3 (66) 
d,(44,/77) COs x+- COS y 
¥(4¢6,/7) 7 Sin o ss 
and it he = ; (67) 
I, (4¢5 Tr) COS ~-+-COS y 
Also we have, from (19), 
L/H 2 sin y/(1-+-cos «cos y). (68) 


By using the above limiting process the values of the lift and moment 
acting on the aerofoil can be deduced from the general values. We find 


Y 7 a: 
* (H?+ 2?)+ H,(H,—2H,)- 
i sity oo Oe 
a7 (Het 3H,(H,—2Hs)}}. (69) 
4 0, ia + 4Y/sin(a+y)/sin y, (70) 
r aH? 
= : H2+-7?), (71 
}oU% — — 2H, Hisinty” 2”? | 
where 
H #4(4¢, 1) (4d, Tr) a(sin y sin x) (72) 
Iy(3)/7) —94(3¢9/77) COS ~-+ COS y . = 
H, — 2H,{21(¢r/7) _ 9 8a(3d1/7)) (73) 
: \3(d, 77) (34, 77)| 
ny — ORBLE 4m) B10 , 
| PZ(O)PHS(, T dy) TT} (0) } 
g/3 Cc ak 
H, Pr (OAs if (75) 
4(4¢,/7) 
H. IF (O)9S(2%s 77) P2(0)_ (76) 
> B(0)F(4d,/7) 7, (0) 
H?(0)5 1)9,(bo/7 
and H, 1 (0), (P1/7)9($9/7) 77) 


FRE ($1 + $2)/7} FH 3($1— $2) 7} 


It can easily be shown that, except for the difference in notation, these 


expressions for the lift and moment are in accord with those of Green. 


12. Numerical discussion 

Green has carried out some numerical calculations for the case of a 
particular circular-are aerofoil with camber approximately equal to 0-097 
in the comparatively simple cases when the centre of the circular arc lies 
on the wall and when the circle on which the are lies touches the wall. 


Two different values (i.e. 0° and 5°) have been taken for the angle of 


incidence of the aerofoil. 
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We have carried out detailed numerical calculations for three different 

















circular-are aerofoils, taking y to be equal to 5°, 12°, and 22° respectively; 

these make the cambers of the aerofoils approximately equal to 0-022, 0-053, 
| and 0-097 respectively. The value of the angle of incidence has been taken 
to be 5° in all cases. Our numerical calculations cover not only the simple 
cases considered by Green} but also the more general cases. The results are 
shown in Tables III, IV. and V, and they are also shown 
Figs. 10 and 11 


In these tables and figures the suffix 0 indicates the lift and moment 


graphically in 


- 


coefficients for a circular-are aerofoil in an unbounded stream. In Figs. 10 


TABLE IT] 





Camber o 0-022 (y 5°), a 5 
2 y, iPaGeuzy| TT, 
5445 I 
47 5416 1947 
Q 53901 *QQO!I 
SOI 532 771 
I j 5410 g940 
6< 5677 1-0426 
I 5545 1°07 35 
5884 7519 1-3809 


Camber o 0-053 (y 2”). x 5 

o (Pidpuy)| TP 
341 I 

| 7209 )O2 

ae 968 

61 482 

41 1456 

711 O94 

} 7512 1-004 


TABLE V 


C‘amber o O-097 | 


Y 2s ) x D 
D/(d; rT, 
I 4 I 
- 
} 
74 ) 
8755 O38 
14 527 15 
The values for the aerofoil wit amber equal to 0-097 approximately in these simple 
s have been recalculated, and we have obtained the same values as Green’s. 


railing edge of the are touches the wall. 
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S. TOMOTIKA, 





and 11 the results for the case of a flat plate shown by dotted lines are also 


drawn for comparison. 
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From these tables and figures it is seen that, for sufficiently small values 
of the camber, as well as for small values of the angle of incidence « 


(e.g. a = 5°) usually adopted in practice. as the circular-are aerofoil 
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approaches the bounding wall both the lift and moment coefficients first 
decrease and then increase to values greater than the values for a circular- 


are aerofoil in an unbounded stream. These results are similar to those 


E - 
i 
J 
; 
/Y=0| 
, (plane aerofoil) va 
yY=5 


mus 9 5 - 
4 “ 3 4 
‘ 
- 
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} ee 
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Kre. 11 
lor a flat plate (4, 5. 6). For greater values of the camber, however, the 
iit and moment coefficients are decreased due to the presence of the wall 
ilues { | ] ae P 
rv all values of / H, even when the aerofoil is placed at small angles of 
A } . ° 
- Incidence as in practice. 
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From Figs. 10 and 11 it is seen further that the effect of the wall upon 
the lift and moment of an aerofoil is greatly modified by its camber. Thus. 
for large distances of the aerofoil from the bounding wall the effect of the 
camber is to decrease both the lift and moment found for a flat plate as 
compared with the lift and moment in an unbounded stream. As the 
aerofoil gets nearer to the wall, both the lift and the moment eventually 
increase, as previously mentioned, above the values for an unbounded 
stream, for small angles of incidence occurring in practice, but the rates 
of increase in the lift and the moment are considerably decreased owing 
to the effect of the camber of the aerofoil. 
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THE RESISTANCE OF A RECTANGULAR METAL PLATE 
WITH AN INTERNAL ELECTRODE 


By 8. D. DAYMOND (Mathematical Institute, The University. Liverpool) 


d 5 September 1950] 


SUMMARY 


\ steady current enters a rectangular metal sheet through a small circular internal 
ectrode and flows out through an electrode which coincides with the rectangular 
boundary The problem of finding the effective resistance is essentially that of 
determining a harmonic function @ which is constant over each electrode. The 
ppropriate complex potential function is expressible in Jacobean elliptic functions ; 
the #-funct . however, are found to be more suitable for the numerical work 
olved, and are used below. Altogether there are four parameters (two of these are 
ordinates of position of the inner electrode) ; numerical values of the resistance are 
ound refore, only for certain special but representative values of these para- 
eter \ brief discussion of the numerical results occurs at the end. 


1. THI nethod of images 


s employed below to determine the effective 
resistance of a uniform rectangular metal plate, when the steady current 
enters through a small circular electrode within the plate, the other elec- 
trode coinciding with the whole perimeter of the rectangle. The solution 
of this problem is a preliminary step in an investigation of induction effects 
nsuch a plate when the current varies with time, and when the electrodes 
are disposed as they are for the steady-current case explained above. Near 
the stationary values of the variable current, where the effect of self- 
inductance is negligible, the conditions are practically those for a steady 
current. The resistance of the plate, which we set out to determine here, 
is then of some interest and importance. This type of problem is also of 
some interest to geophysicists who are concerned with the distribution of 
potential when an electrode is pushed into the ground. 
In some actual experiments which were carried out on induction and 
discharge effects, the current entered by a small electrode at the centre 
{the plate and left by a very large number of conducting wires attached 
to the edge of the plate and distributed over the whole of its perimeter. 
For purposes of calculation this can be assumed to mean that the current 
Hows out perpendicularly through the perimeter, i.e. that the whole of the 
perimeter is an electrode. This boundary condition is therefore assumed 
in the following determination of the resistance of the plate to the flow of 


steady current. 


[Quart. Journ. Mech. and Applied Math., Vol. IV, Pt. 1 (1951)] 
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Notation 

a,b (b >a) the dimensions of the rectangle. 

7+ = ib/a_ the modulus of the #-functions. 

q exp(inr exp( mb/a). 

t,o the constant thickness and constant specific conductivity (respec 
tively) of the plate. 

x.y Cartesian coordinates with origin at the centre of the plate and 
axes parallel to the edges (Fig. 1). 


ty. Yq the coordinates of the centre of the inner small circular electrode. 


L+-1Y, 2s Lot lo; 
v E+in Z/20, Vo Eo tino Z| 20. 
w—=go+is the complex potential function (the current density is pro 


portional to the gradient of ¢). 
» the potential of the inner electrode. 
ce the radius of the inner electrode. 
) c/2a. 


R the resistance of the plate. 


In section 2 the function w, and subsequently ¢, and the expression 


2rotR, are derived for such a system in which both the plate dimensions 
and the position of the inner electrode are arbitrary. The expression for 
dy iS & power series in 4, and terms of the order 6° and smaller are ultimately 
neglected. 

| An expression for w can be obtained, by a conformal transformation, in 
terms of Jacobean elliptic functions: it is, in a usual notation, 


w 2 In{(t—t,)/(t—t,)}. 
where ¢t = sn{r.K(2z+-a)/b}, 
and K'/K 2a/b determines the modulus /. For carrying out the numeri 


cal work, however, this form is less convenient than that derived below. | 
2. The current density in the plate is i ogradd, and, since in the 

space between the electrodes divi is zero, therefore 
V2dh 0. (1) 


The total current flowing out of the plate is 
t | (i-ds), 


where the integration is taken round the rectangle. The integral is equal to 


f 
? ds at| as], 
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by a property of conjugate functions. The square bracket indicates the 
nerement which results from one complete description of the rectangular 
oundary. Taking the potential on the latter boundary to be zero, the 


value ot the effective resistance is 
at)! dy [yp]. (2) 


The function w. whose real part must satisfy (1) and the simple boundary 
mditions, is now determined by using the analogy between this problem 
nd that of the two-dimensional conductor at zero potential under the 
nfluence of a unit positive line charge passing through the interior position 


) occupied by the electrode (see Fig. 1). 
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Ww 
In this electrostatic problem the given charge and the image system 
the nsist of line charges 
| through the points v Vo tm-+nt 
l i tg “ 2 
l 7 nn {o) 
| through the points 1 Vo ttert+tm-+nr 
| Vo } m-+-NnT, (4) 
ere vy is the complex conjugate of vy. and m. n 0, +1, +2..... 
Che required function w is —21n{f(v)} where f(v) has simple zeros at the 
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points (3) and simple poles at the points (4). Now 3,(v 7) has simple zeros 
at the points v = m-+-nr (see 1, p. 255). Hence 

F(v—v_)P,(v + V9 } Lr 


) 
F, (v—vy—4r)3,(v-+0,—4) 


» 


cu In 


3, (v—v_)P3(v+Vp) 


In- a : 
a,(v Vo)Po(v +-Vo) 


B,(v—vp)Pa(2v9-+v—vp9) 


I(t 2N9+v—Vp)Po(2Eq+-v Vo) 


In 


With the help of (1), p. 252, equations (1) to (4), and Taylor’s theorem 
we now express the right-hand side as a certain function of v—v, and a 
series expansion in v—v,. In the series the coefficients of v—v, and (v—v,)* 
are simplified by the use of (2), p. 100, equations (4). Also, since the third 


and higher powers of 6 = c/2a are being neglected, the potential ¢, on the 


inner electrode is determined by substituting either 7 No: € = &y+8. or 
e= §,, 4 Ho +5, and retaining the real part of w. Thus 


1 Hd. cos 27€, 
dy in( : ) In J 


7D, sin 7d 
2> Q,,|¢'"(1 — cos 2mzd) + (— 1)"(1 —cos 4mz€, cosh 4779) 
I 
(cosh 4ma7n,— 1)+(— 1)™q""(1 — cos 4mz,) | 
tm 


+75| 1 tan 27€, > (—1)"-'mQ,, sin 4mzE,(q’" — cosh tmzn9)| 
i 


47°5?| 1 sec 2n€, 
oO 
> mQ),, {cosh 4many+(— 1)” cos 4m7é,(q"" cosh tman)}| 1 ()(§%) 
(5) 
where # represents #(0/7) and Q,, = q'"/m(1—q?"). 
The function y for a line-charge of unit strength, through the point z = 2, 
is —2arg(z—z,). In this case lh is equal to —4z or zero according as the 


contour encloses z, or does not. When, therefore, the contour is the rect- 
angle itself, only the complex potential arising from the given line-charge 
will be effective in the determination of |] as the images are all outside 
the rectangle. Thus [¢] -47. From (2) the resistance of the plate is 
(2zrat)-1 345. 


3. The largest value of q is exp(—7z) = 0-04321... when b/a is unity, 


and qg decreases from this value to zero as b/a increases from unity to 
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infinity. The function In(?\i, 7%, 3,) also decreases steadily to zero as b/a 


increases in this way: some pairs of corresponding values are: 


ba 1-0 1-5 2-0 2°5 3 

In(, d, 7H, 7,) 0165804 0-035614  0-007456 = 0-001552 — 0-000323 

\ll of the series in (5) tend to zero as b/a tends to infinity. The rate of 
convergence of at least one of them is, however, slow, particularly in the 
case of a system where the plate is square and the inner electrode is fairly 
distant from the centre 


In the first table below are some values of 2zotR. for the range 


| bh a x 


ind for the positions of the inner electrode represented by 





; Y¥. 2a = 0 10€, = 20a,/a = 0, 1, 2, 3, 4, and 5. 
3 } 5 
+75 F755! 7124 $°0397 
5 6500 3*O055 3°5 320 
45799 | 4°5174 
53 777 37119 
1437 | 49053 | 4°8397 
j *OO2I 37339 
{ #9100 | 4°5443 
{ 8074 | 3°7384 637! 
}QO115 $5453 7473 
494 5055 7394 0363 
J 5 HOI4 $O49E 4°75 2K 
158 3°7 506 54 





rhe two values of 270tR occurring under each set of values of b/a. a», and 
j). are such that the upper and lower values correspond respectively to the 
alues 1/450 and 1 150 of 5. The latter have been chosen partly because 
they are of the right order of magnitude and partly as a convenient pair 
of values (in the computational sense) to serve as a basis for the com- 
parison of two different sets of numerical results. 


When x, and y, have fixed values the convergence of each of the series 


( 


n (5) is slowest when b/a |. In this case, because of the symmetry of the 
square plate about a diagonal, the convergence of each series must be 
slowest when x) = Yo. (This is also confirmed by some numerical results 


not included here.) 


The short table below gives pairs of values of 27otR for the case of a 
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square plate when the inner electrode has different positions represented by 





A = 202,/a = 20y,/a 1, 2, 3, 4, and 5. 
A I 2 3 rl 5 
27rat hk 1°7808 4°7 284 pO412 4°5177 | 4°3539 

"6503 | 3°6267 375350 | 3°4129 | 3°2472 


For any value of b/a the resistance decreases as the inner electrode is 
moved from the centre towards an edge. It decreases also in the case of 
the square plate when the inner electrode is moved from the centre towards 
a corner. The effect on the resistance is relatively greater, however, when 


the radius of the inner electrode is increased; it is true, for all values of } a 


and for all the electrode positions considered here, that the resistance is rot 
reduced by about one-quarter when the small radius is trebled. On the - 
other hand, a movement of the electrode from the centre of the plate to pe 
a point half-way between the centre and one of the longer edges (a move- on 
ment of 30 diameters when 6 1/450) causes a fall in the resistance of less bo 
than one-twelfth. me 
se] 
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rds 
hen SUMMARY 
h his note describes one nd two-parameter families of solutions of steady 
By des tationally-symmetric viscous flow. The solutions are such that the Navier-Stokes 
equations reduce to ordinary differential equations in a single position variable. 
- The one-parameter family represents flow which is rigid-body rotation at infinity 
ind over a plane through the origin; the solution given by von Karman in 1921 is 
ve me member of this family. The two-parameter family represents flow which is rigid- 
eas body rotation over each of two planes at a finite distance apart. The case of large 
Reynolds number is particularly interesting, since the two bounding planes are then 
separated by a region of rigid-body rotation and translation in which viscous effects 
are negligible. 
Introduction 


T. v. KARMAN (1) has pointed out a simple solution of the Navier-Stokes 
equations of motion which describes the steady flow of a viscous fluid in 
a semi-infinite region bounded by an infinite rotating disk. This ‘solution’ 
is not yet given analytically, since one is left with two ordinary non-linear 
differential equations in a single independent variable which must be solved 
numerically, but to have carried a solution of the Navier-Stokes equations 
even so far by exact analysis was (and still is) something of a novelty. 
Moreover, the solution has the very interesting property that it is also a 
solution of the appropriate boundary-layer equations, the terms neglected 
in boundary-layer theory being identically zero for this type of motion. 
The purpose of this note is to show that there are one- and two-parameter 
families of solutions having the particular mathematical simplicity of Kar- 
man’s solution; Karman’s solution is one member of the one-parameter 
family. In Kaérman’s problem the flow far from the disk is assumed to be 
wholly normal to the disk and to be induced by the rotation of the disk. 
In general, if other conditions far from the disk are assumed, the particular 
simplicity of Karman’s solution is lost, but it will be shown below that the 
simple form of the solution is retained if the fluid at infinity has an arbi- 
trary uniform angular velocity y about the axis of rotation of the disk. If w 
is the angular velocity of the disk, there is found to be a solution for each 
lue of y/w between © and *. Karman’s solution corresponds to 


[Quart. Journ. Mech. and Applied Math., Vol. IV, Pt. 1 (1951)] 
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y/w = 0, whereas y/w Foo describes a flow which is rotating uniformly 
at infinity and which is bounded by a stationary disk. In this latter case 
there is thus the opportunity of describing quantitatively the tendency for 
particles of sugar to migrate to the centre of a cup of tea which has been 
stirred. The qualitative explanation (2) of this phenomenon is well known, 
of course; the new point is that the flow away from the disk induced by the 
rotation is uniform over planes parallel to the disk, just as the flow towards 
the disk is uniform over these planes in Karman’s problem. 

A two-parameter family of solutions of the same simple type describes 
the flow between two parallel infinite disks which are rotating about the 
same axis with different angular velocities; in addition to the ratio of the 
angular velocities of the disks, the Reynolds number based on the distance 
between the disks is now a relevant parameter. A numerical method of 
determining the flow field in the special case in which one disk is stationary 
and the Reynolds number of rotation of the other disk is small has recently 
been described by Casal (3). 

A family of solutions of related type has been described by Miss Hannah 
(4). This family is obtained by combining the flow towards the disk pro 
duced by a source at infinity on the axis of rotation with the rotating flow 
induced by the disk. Karman’s solution is obtained when the source-flow 
is made zero, and, at the other limit, if the disk is stationary the viscous 
stagnation-point solution described by Homann (5) is recovered. This 
family of solutions will not be included in the following discussion. 


The governing equations 


If v,, vg, v. are velocity components in the directions of increase of cylin- 


drical polar coordinates 7, 6, z, and p is the pressure, the Navier-Stokes 
equations of steady motion of a fluid of density p which is symmetrical 
about the axis r = 0 are as follows: 
ov ov. vz Le v,,. 
. ) ‘ Fs 
v. — +9, — e xe. Poy V2v,——"], (1) 
cr " C2 r por ” 
CVp Ovp U,V” a Vp 
v, OQ, y 9 “r'd 1 V2V9- - :) (2) 
“ ni -~ . 2 
C7 Cc } 
ov, ov, l Op ; F 
v, —+v, — u 1 pV2v., (3) 
cor " OB p Oz c 
oe l € oC 
where V2 = hiow 3 , 
‘ 9 ‘ 9 
Or? ror 02 
The continuity equation is 
l orv, , dv, 
— = — 0 (4) 
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The plane z = 0 is identified with the plane of a uniformly rotating disk 
(angular velocity w) so that one set of boundary conditions is 

v v 0, ve=—or, at z= 0. (5) 

The other set of boundary conditions will depend on the problem under 

discussion; if the fluid is unbounded in the z-direction the conditions to 


be assumed are 


v1, > 0, Va>y,7, @ 270, (6) 
whereas if there is a second rotating disk at z = d the conditions are 
v v. 0. V9 Yor, at Zz d, (7) 


where y;, vy, and d are disposable constants. This excludes Miss Hannah’s 
family of solutions, which give a radial velocity at z x, but it includes 
the one- and the two-parameter families mentioned in the introduction. 

The simple property of Karman’s solution and of the solutions sought 
herein is that the flow normal to the disk is uniform over planes parallel 
to the disk. That is, 


v v.(2), 
and as a consequence, from (4), and assuming that v, is finite at 7 0, 
r dv. 
v, ~ (3) 
2 dz 
Equation (3) can then be integrated to give 
dv, , 
Pi y—*— }y*+ I(r). (9) 
dz es 
From equation (1) we find that the arbitrary function II(7) satisfies 
ldil(r) 2% . ; . 
“ — function of z only. 
r ad y2 


The boundary condition (5) then requires 

I(r) = $7?(w?+-c), (10) 
where c is a constant, so that 
Ug)? function of z only. (11) 


t 


Equations (1) and (2) then become 


l dv.\? d*v Vp\" v dy. 
Ly 0 = vp 7-~C = 2 
(5 ae) = dz ("") ee 2 dz (12) 
and dv. UG v d(v9/r) v d*(vp id (13) 
dz 1 dz dz 


In the case of the semi-infinite fluid with boundary conditions (5) and 


6), the asymptotic form (as z -> 00) of (12) is 


yi—w? (14) 
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and only five boundary conditions are needed for the solution of (12) and 
(13). These are supplied by (5) and (6), the uniform axial velocity V at 
2 = © being determined as part of the solution. In the case of the two 
rotating disks all six of the boundary conditions (5) and (7) are needed to 
solve (12) and (13) and to determine the constant c. 

The form of solution which has been assumed is still valid when there 
is a uniform suction through the surface of either or both of the disks, the 
boundary condition v, = 0 being then replaced by v, = a prescribed con 
stant. Suction through the surface of the disks will give rise to some 
interesting modifications of the flow patterns described below, but will not 
be considered further. 

The variables may be made non-dimensional by using (vw)! as a reference 
velocity and (v/w)! as a reference length, the direction of positive rotation 


being so chosen that w is always positive. Put 


r (v/w)*n, Z (v/w)*Z, 9 (vw)? ng(C), v. (vw)*h(Z), 


in which case equations (12) and (13) become 


Lh’2—Lhh" —¢? ~ i ‘ hh”, (15) 

w* 
gh'+-q’h q’, (16) 
where dashes denote differentiation with respect to ¢. The boundary con 
ditions are now h h! 0. g a oe 0. (17) 

and, for the semi-infinite fluid (in which case ¢ = yj—w?), 

h’+0, g>y,/w, as (>, (18) 

or, for the fluid between two rotating disks, 
t= = §, g= Yew, at C (=. (19) 

> 


Numerical integration of equations (15) and (16) by the method used 
by Miss Hannah (4), and by Cochran (6) in the case of Karman’s problem, 
would probably be feasible, but tedious. Interest lies more in the general 
form of the solutions than in the numerical details so that only general 
observations about the streamlines in typical cases are presented here. 
The author has no evidence for the conclusions stated, other than that 


mentioned explicitly. 


Streamlines of the flow bounded by a single disk 
Consider first the family of solutions given by the boundary conditions 
(18). There is one member of this family for each value of y,,w between 


© and +, and for each member there is an appropriate value of J’, the 
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q| axial velocity far from the disk. The members of the family may be divided 
- into three classes within each of which the streamlines have much the same 
"0 appearance. 
"to Class (a) t Y)/w | 
| The rotational velocity at the disk is, in this case, smaller than at any 
a ther point in the field so that the inward radial pressure gradient imposed 
he 
n 
ne 
ot 
2b 
ce 
on 
dD 
16 
yn 
7 < 
0 
disc r 72 
Ss 
| if 
- oo 6S Vw) > | 
Fic. 1 
sed . . > . . . . . 
by the fluid at infinity is more than sufficient to keep the fluid near the 
Pm ] . . . . . . . 
, lisk moving in circles. Hence there is a radial flow inwards at points near 
‘Ta 1 . ° . : ~o . ° 
. the disk and an axial flow away from the disk (see Fig. 1; in this and other 
ral . . . . . . 
| figures the streamlines refer, of course, to components of the motion in an 
ere " " > 
: ; axial plane only). At one end of the range y,/w = 0, corresponding to 
lat : ; ee : : 
’ stationary disk and giving an approximation to the tea-cup flow, and 
it the other extreme y,/w = 1, corresponding to uniform rotation of the 
whole fluid with no axial motion. 
ions | Class (b). 1 > y,/w > 0 
een | The angular velocity at the disk is here a maximum and the rotation is 
the everywhere in the same direction, so that the axial velocity is towards the 
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disk (Fig. 2). The disk acts as a centrifugal fan, throwing fluid out radially 
and drawing it in axially. The extreme case y,/w = 0is Karman’s problem. 
Note that the streamlines of members of this class are not obtained simply 
by reversing the streamlines of members of class (a). 

Class (c). 0 > y,/HW > —a 

The disk, and the fluid far from the disk, rotate in opposite directions, 
in this case so that at some value of € (i.e. of z), g(f) = 0 and the fluid 


x, 














there has zero angular velocity about the axis r = 0. The angular velocity 
of the disk is thus greater in absolute magnitude than that of the fluid near 
it, and we anticipate that the radial velocity will be outward in the neigh- 
bourhood of the disk (i.e. the disk acts locally as a centrifugal fan) but 
inward elsewhere. This conclusion is supported by the following rough 
investigation of the form of the functions g(Z) and A(Z). At € = 0 we have 
g = | and, it may be assumed, g’ < 0. When Z is large, g is asymptotic to 
the constant values y,/w so that the function probably has the form shown 
in Fig. 3. Now equation (16) can be written in the form 


4 


hg=—[Gat or =— [a+ ok, (20) 


J 9 (vw)! 


” ie a 




















o> Vi fw > - eo 
Fic. 3 






























































36 G. K. BATCHELOR 


whichever range avoids the singularity at the point where g = 0. On the 
supposition that g” > 0 for all ¢, as in the sketch in Fig. 3, the variation 
of h required by these equations is sketched in Fig. 3 and the streamlines 
are shown in Fig. 4. Equation (16) shows that the value of ¢ for which 
h = 0 must be at least as great as that for which g = 0. Apparently the 
critical plane on which the axial velocity 4 vanishes, which may not 
coincide with the plane of zero angular velocity, divides the flow field into 
two self-contained regions. The dividing plane coincides with the disk 




















when y,/w = —oo and moves away to infinity as y,/w increases from — 
to U0. 
disc 
a 
+ 

a) 

_ 4 

—_ . 

disc 


co > yz,/w > 0 
Fia. 5 

Streamlines of the flow between two rotating disks 

The streamlines in typical cases can again be sketched from elementary 
considerations. There is one solution for each value of the parameter y./w 
between —oo and +o and for each value of the parameter d?w/v between 
0 and +o. Variation of the parameter d?w/v has the effect of varying the 
extent of the region of rapid change of angular velocity, i.e. of controlling 
the boundary-layer character of the flow. The form of the streamlines does 
not vary radically with d*w/v, so that the whole two-parameter family may 
be divided into two classes in which y,/w takes opposite signs. 

Class (a). +0 > y,/w > 0 

In this case the disks are rotating in the same direction and the magni- 
tude of the angular velocity of the fluid varies monotonically with ¢. The 
radial velocity will be inwards near the slower rotating disk and outwards 
near the faster, which acts as a centrifugal fan, and the streamlines will 
be approximately as in Fig. 5. The extreme case y./w = 0 has been 
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considered recently by Casal (3), who showed how the equations (15) and 
(16) can be integrated by expanding the velocities as power series in 
d?w/v which are convergent when d?w/v < 0°17. 

An interesting situation arises when the Reynolds number d?w/v becomes 
very large. The effect of the no-slip condition is then confined to thin 
layers near each disk and the flow outside these layers is approximately 
as for a frictionless fluid, i.e. the angular velocity in this region is approxi- 


mately independent of ¢. The streamlines are sketched in Fig. 6. The 
disc 


oe am as 
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co > 2 /” >0,d’w/y —>oo 
Fic. 6 


uniform angular velocity in the interior of the fluid will presumably have 
a value such that the axial flow away from the slower rotating disk (con- 
sidered as a disk rotating in a semi-infinite fluid which has constant angular 
velocity far from the disk) is just equal to the axial flow towards the faster 
rotating disk (considered in the same way). Thus if V(y,/w) is the axial 
velocity far from a disk which has angular velocity w in a semi-infinite fluid 
which has angular velocity y, far from the disk, the angular velocity Q in 
the region between the two disks at large Reynolds number is given by 

V (QQ/w) V(Q/y.). (21) 
The consideration of a single disk in a semi-infinite fluid showed that there 
will be a solution of this equation provided that, as already supposed, 0 
is intermediate between w and yp. 


Class (hb). O Vo/W x 
The disks now rotate in opposite directions and there is some plane 
between the disks on which the value of vg is zero. Thus the radial flow 


in the neighbourhood of each disk will be outward, with an inward radial 
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flow in the interior of the fluid, as sketched in Fig. 7. As in the case of the gene’ 
semi-infinite fluid there will be a division of the flow into. two self-contained Whe 
regions, the dividing plane on which v, 0 being not necessarily identical vani 


(so far as can be seen without detailed numerical work) with the plane on angu 
which vg = 0. The dividing plane coincides with the upper disk (angular dete 
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disc ¥2) 
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>... 
disc (w) 

} 
O> y,/w > -— © pas 
Fic. 7 tuit 
tha’ 
velocity y,) when y,/w = 0 and moves down to the lower disk as y,/w and 

decreases from 0 to —oo. 

The case of very large Reynolds number is again very interesting. The 7 
effect of viscosity falls off rapidly with distance from each disk and outside not 
thin layers near each disk the angular velocity is uniform and the stream- disl 
lines in a plane through the axis are axial. The uniform axial motion out- vel 
side the boundary layers will normally be from the slower to the faster | the 

y tay : } 
rotating disk. However, the magnitude of the angular velocity decreases pro 
with distance from each disk and each disk must therefore act, locally at | — son 
least, as a centrifugal fan. In the immediate neighbourhood of each disk cha 
the axial velocity must therefore be towards the disks, and the axial cha 
velocity must change sign somewhere between the disks. The means where- spe 


by this can happen have already been explored; the flow in the boundary ) tw 
layer near one disk (the slower) is presumably as sketched in Fig. 4, and | _ the 
near the other (the faster) as in Fig. 2. The streamlines thus have the sta 
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general shape shown in Fig. 8, which is drawn for the case |y.| < |w!. 
When ‘y,| decreases to zero the outward radial flow near the upper disk 


vanishes and the flow reverts to that obtained by inverting Fig. 6. The 
angular velocity Q in the region outside the two boundary layers is again 
determined by the equation 


V (Q2/w) V(Q 


5 
Ys 


> 


), 


where V has the same meaning as in (21). Since Q has the same sign 


dise (22) 























y Y y ang. vel. 2 























disc (w) 
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O> lu >- mo, A o/y —> © 


} Fic. 8 





is whichever of w and y, has the greater magnitude (so one believes in- 
tuitively), consistency of this equation for Q with Figs. 1, 2, and 4 requires 


that : 
ass if Vo/W l, then Q « Q)/w 3 


and reciprocally, 
if \yo/w l, then 0< Q/y, <1. 

There does not appear to be any mathematical reason why there should 
not exist a solution for which the boundary layer on the faster rotating 
disk is the one in which the angular velocity changes direction. The axial 
velocity in the region outside the two boundary layers is in this case from 
tl 


probably be unstable, unlike that described above. The situation bears 


e faster to the slower rotating disk. However, the solution would very 


some resemblance to two-dimensional viscous flow through a diverging 
channel, and the analogy is closer if we imagine the two plane walls of the 
channel to be moving in their planes in the direction of flow at different 
speeds. There will be a region of reversed flow near one of the walls and 
two different solutions will be possible. But again it is probable that only 
the solution which gives reversed flow near the slower moving plane is 


stable, 
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The case y, w is singular, because there must then be a possible 
distribution of velocities symmetrical about the mid-plane (just as there 
is a singular symmetrical solution for flow in the diverging channel when 
the two planes have equal speeds—or in particular are stationary). Sym- 
metry about the mid-plane implies that the boundary layers on each disk 


are mirror images of each other and that the axial velocity immediately 
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outside each boundary layer is in each case towards the disk. Somewhere 
in the interior of the fluid the axial velocity must change appreciably and 
must reverse its direction. This‘cannot happen in the absence of a strong 
viscous effect, so that there is apparently a transition layer of rapid change 
as sketched in Fig. 9. The axial velocity and the angular velocity both 
change sign within this layer, while outside it (and outside the disk 
boundary layers) they are constant. 

This singular solution may not be realizable experimentally, of course, 
but it has some intrinsic interest. For instance, it indicates that there is 
yet another solution of the type considered herein. The central transition 
layer is not directly connected with the two boundary layers (and in the 
limit of infinitely large Reynolds number, its position is arbitrary, provided 
it does not overlap with either of two of the boundary layers) and can be 
regarded as a transition region between two semi-infinite masses of fluid 
rotating with equal and opposite angular velocities =Q. Such a flow is 
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described by equations (15) and (16) (with ¢c = Q?—w?; note also that the 
reference angular velocity w ought now to be replaced by Q) with the new 
boundary conditions 





h=h" = 0, g=0, at (€=0, 
and h’>0, g>FQ/w, as €> Fo. 
, a disc (-w) 
- = 
Boundary 
layers 














dise (w) wo ol 


Distribution of _anguiar velocity between the dises. 
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If this solution gives an asymptotic axial flow U(Q/w) away from the 
transition layer, the condition which determines the angular velocities -Q 
in the non-viscous regions of the flow in Fig. 9 is 
U(Q/w) V(Q/w), 
which will have a solution provided Q) < |w!. The angular velocity of the 


fluid between two disks ev idently varies with z in the manner shown in 
Fig, 10. 
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THE VELOCITY DISTRIBUTION IN THE LAMINAR 
BOUNDARY LAYER BETWEEN PARALLEL STREAMS 


By R. C. LOCK (Gonville and Caius College, Cambridge) 


[Received 13 June 1950] 


SUMMARY 


A method is given for obtaining the solution of the laminar boundary layer 
equations for the steady flow of a stream of viscous incompressible fluid over a 
parallel stream of different density and viscosity. An approximate solution is also 
obtained by means of the momentum equation. It is shown that the solutions depend 
only on the ratio U,/U, of the velocities of the two streams and on the product py 
of the corresponding density and viscosity ratios. Numerical results are given, in 
the case where the lower fluid is at rest, for four values of py, and also when py = 1, 


for one non-zero value of the velocity ratio. 


1. Introduction 


[IN a paper on the stability of the flow of a stream of fluid over a layer of 


the same fluid at rest, Lessen (1) has obtained the velocity distribution 
of steady motion in the free laminar boundary layer separating the two 
streams, using a method equivalent to that of Blasius for the boundary 
layer on a flat plate. In the present paper the more general problem, when 
the two fluids are of different densities and viscosities, is considered. 
Sir Geoffrey Taylor, in an unpublished note, has already given a simple 
approximate solution by von Karman’s momentum integral method, but 
it was thought worth while to investigate the accurate solutions of the 
boundary layer equations for this case, since it is hoped later to consider 
the problem of the stability of the same motion. 


2. The boundary layer equations 


We consider the two-dimensional motion of a stream of fluid with 
velocity U;, density p,, and viscosity ,, over a parallel stream with velocity 
U,, density py, and viscosity p,. Both fluids are assumed to be incompres- 
sible. 

We take the axis of x to be horizontal, in the direction of motion of the 
free streams, and the axis of y to be vertically upwards. The origin is taken 
as the point at which the two fluids are supposed first to come into contact 
(see Fig. 1). The corresponding components of velocity are wu and v. 


[Quart. Journ. Mech. and Applied Math., Vol. IV, Pt. 1 (1951)] 
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place in a layer of small thickness, and that v is everywhere small compared 
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with uw, the boundary layer equations are 
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u — 0 Vy = (1) 
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it10n 
two for the lower fluid, where v, and v, are the kinematic viscosities of the two 
= fluids. 
ee The equation of continuity is 
red Cu cv ‘ 
0, (3) 
nple Ca CY 
but so that there exists a stream function w& such that 
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cy Cx 
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Equation (1) then reduces to The 
- eae and tl 
2 afi Lf, a“f, (0. (9) | Wem 

dy? : dn? 


In the lower fluid it is convenient to use a different variable », defined which 





by U7. 1 The 
Ne : J; Y, (10) 14) w 
Vg t 
that (" : 
so é No nn. 
‘ = " Thi 
and to put us (vy U, x)*fs(79). (11) For, 1 
ga since 
Chen : U, fin), (12) | to the 
I 12 4 ri or ) : 
v= = (ne fa(n2)—falna)} (13) 
ee ee 
and ee v,( fel) (14) In | 
cy Vot 


Equation (2) then reduces to 


3 d2f 
9 d*f, P If, (0). (15) 


2 —“= 4+ f, —$ 
dy dn Usi 
Boundary conditions 
The boundary conditions at infinity are ' 
Irom 
u>U, as y,>+0 
and u>U, as n> —, 
so that ti ~ | as Ut a t+- OO (16) No’ 
| 
~ i . 
and h->— as n> —O. (17) 
U, , 
Since the motion is steady, the interface between the fluids is the | ‘mila 
streamline % = 0 which passes through the origin, and is therefore given f(y?) 
MN 
by ; 
: h=h= 0. | The 
If »$ is the value of y, such that f.(y$) = 0, then we must have also The 
Silat) = 9, 
i 
Vo\5 
where nf = | ia n> Iti 
vy conti1 


(It will appear later that »? = 72 = 0 unless p, = py.) interf 
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The other conditions to be satisfied at the interface are that the velocity 
nd the normal and transverse components of stress should be continuous. 





( We must therefore have eo - ' 

’ fila) fo (3), (18) 
fined | which ensures that both uw and v shall be continuous, even if 7}? + 0. 
The tangential stress is u(Cu/Cy) approximately; using equations (8) and 
(10 14) we see that this will be continuous provided that 
, } an 

Pi ¥i Si (71) p23 fo (72). (19) 
This result may also be obtained from considerations of momentum. 
1 For, if the boundary between the two fluids is given by y = y,(x), then, 


ince there are no solid boundaries and therefore no frictional resistance 
(12) | tothe motion as a whole, 
iyo 
3 |} p,u(U,—u) dy+ | pou(U,—u) dy 0. 


x 


(14 In terms of the variables 7, and 7, this becomes 
Ne 
— Pg [Us ge 
Pi) pifi(l—fi) dn 4 Pots |+> —f2) dn = 9. 
‘ 1 
(15 , . 


Using the fact that 


[f2dq=Sf'—[ ff ay 


rom (9) or (15), we get 


Bh U) 
PVilli—hd — Ah | pat = fo—Iefe —2fs - 0. 
l 1 —o“ 
6 Now f; («) 1, fi (@) 0. and it can be shown that 
7 | lim f,(i—f,) 0; 
‘. —_ y sia : 
; the | Smilarly fo (—o) 7 Ja ©) = 0, and lim f,|-*—/f,) = 90, and also 
. l es Uy 
o1ven n’) fol 72) 0. 
Therefore, p, vi f7(7?) p2v3 fs (m2) as before. 
| ow 
| The nor mal stress 18 f 
cov 
Darn pt+2p : 
cy 


It is easy to see that if u(@u/ey) is continuous, then p(év/éy) will also be 
mtinuous. The pressure p must, therefore, also be continuous at the 


nterface. 
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Now the variation in p/(p, VU?) across the boundary layer is known to be 
of order (3/x)*, where 6 is the boundary layer thickness. On the other hand, | 
the difference between the hydrostatic pressures at the two sides of the 
interface is given by 


Pi—Pe2 = ( = Ps) JY 
pil i Pi) l i 


7 ( Pg (2) nt 
Py j U, Uy 


1 is of order unity, this is of order 8/x. 


and, if p, ~ p, and 7 Thus when 
Py < pz it is evident that 7} must be small, of order 5/x, and may be taken 
as a sufficient approximation to be zero, so that the effect of gravity will 
be to maintain the interface horizontal. 

In this case, therefore, the solution of equations (9) and (15) is com- 


pletely determined by the boundary conditions 


fi(~) = 1, (20) 
a U, . 
fa ( 20) U,’ (21) 
fi(9) = f2(0) = 9, (22) 
f1(0) = fo (0), (23) 
and Pi V4 fi (0) = pov} fe (0). (24) 


When p, = pz, however, gravity has no effect and 7 is left apparently 


arbitrary (see 1). This does not in fact materially affect the solution; there 


are certainly an infinity of solutions of the equation 


ff'+2f" =0 


which satisfy the boundary conditions 


but if f and g are any two such solutions it is easy to see that they must 
be connected by a relation of the form 


f(n) 


where b is a constant, so that also f’(7) 


y(n T b), 


g' (n+). 
The difference between the two solutions is, therefore, simply equivalent 
to a shift of the velocity distribution as a whole in the y-direction. 

In a given physical case, however, there must be a unique solution and 
a definite value of 7}. In order to fix this it will be necessary to specify a 
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further condition at infinity. Now it can be shown that the displacements 
of the streamlines at great distances from the interface are 
oF " “) lim (,—f;) 
1 _ 
for the upper stream and 
8% Pe } lim (7 —fo) 
i 

for the lower stream. These displacements will depend on the value of 7%, 
so that one of them may be chosen arbitrarily without interfering with the 
flow in the boundary layer, and the physical conditions will then be fixed 
completely. (If U, = 0 only df has any meaning and this may be chosen 
arbitrarily.) It should be noted that the situation is quite different from 
the case of boundary layer flow along a flat plate, when the presence of a 
constraint at infinity, such as a wind-tunnel wall, will cause definite 
interference 
3. Approximate solutions by the momentum equation method 

sefore discussing the accurate numerical solutions of equations (9) and 
(15) it seems of interest to obtain some idea of the general nature of the 
velocity distributions by means of von Karman’s momentum equation 
2, p. 131). 


The momentum equations in this case reduce to 


8 
70 ; -| . OS 95 
vs (U,—u)u dy, (25) 
Py CY] y-+0 CX . 
0 
for the upper fluid, and 
82 
T Cu c ; - 
[= Vy | (U,—u)u dy, (26) 
P2 “ley -0 Ox 7 
0 


for the lower fluid, where 7, is the skin friction and 5, and 8, are the 
boundary-layer thicknesses on either side of the interface. 
In the upper fluid, we suppose that 


ul 
ok 
ly $,(77), 
1 
' ’ 
where Ni y 
mr . * 
Then equation (25) becomes 
— 
V dé % 
Ly 2) dn* 
: $,(0) 0; ; (d, $5) dy}, 
U; dx . 
0 
1 
» > fe 
2y.% » as , 
whence 5 $,(0) oF | (4, —¢2) dnf. (27) 
l ‘ 
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In the lower fluid, we suppose similarly that 


u 
> = $2(73), 
U; 
where n> _¥ 
Ps O»5 
and then equation (26) becomes 
1 
2,2 , Q9 f 2 * 
77 $20) = 35 | (Ab, —3) dy¥, 
d 
U, 
where A= —. 
U; 
; 
If we write TI, = | (¢$,—¢3) anf, 
0 
2 
I = | (Ad, ~$3) dng, 
0 


Pi yy et Y 
=> = 0) ~ — ° 0 " 
0 5, Pl 5, 2 ) 
we get pol,+f, = 0, 
—— 2 
where 6 = — and p = Pa. 
1 Pi 


the conditions 


(1) = 1, $,(1) = $3(1) =... = 0, 
$3(0) = £7 (0) = ... = 0, 
$,(0) = $,(0), 


By ‘(0 oy cane "(0 
5, Pl ) 5. Pal ), 


$3(0) = $3(0) = ... = 0, 
$o(1) = A, $4(1) = 69(1) = ... = 0 


found from equation (29). 


possibility of using the simple expressions 
$, = 14(c—l)e-7i 
and dg, = A+(c—A)e-™ (with A= 0). 








and eliminate x from equations (27) and (28), making use of the relations 


(29) 


The functions ¢, and ¢, may be chosen so that they satisfy as many of 


as is convenient. The integrals J, and J, can then be evaluated, and 6 


3.1 In suggesting this problem, Sir Geoffrey Taylor pointed out the 
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The boundary conditions at 7f = 1, ny = 1 are now replaced by similar 
conditions at infinity, all of which are satisfied. The conditions 


b,(0) = (0) = c, 


and 54;(0) id,(0), where pu to/,, are also satisfied, provided that 
1 LPe I Me/ Fy i 
¢ a OTA (32) 
o- es 
Then IB | (6,—¢%) dnf 1(1—c?), 
uit I, | (Ads $3) dns y(r? c*). 


Substituting these expressions in (29) and using the value of ¢ given by 


») 
) 


2) we obtain eventually the equation 


pl 25+ u(1+-A)| = pd*[6(1-+-A)4 2pd]. (33) 
This may be written in the form 
6\? l s 2(¢ T l 1)! 
| ) (= (0/2) d' (34) 
bu py | (oO p)(1 t A) | 2A] 


which is suitable for solution by successive approximations, the first 


spproximation when py is large being 


( I+A\i ., 

: }° if A 0 
Le 2A pL, 
re) 


pL) if A 0. 


[L 
i 


The velocity uw, at the interface is given by (32); 


' Wo, A+ Olu) 
1S 67 3 
U; 1 +-(0/p) 
l . U] O/[L) . co ° 
so that up, : \O/t U,+ U\(d/p) if pp is large. 
l (0/p.) . 


The skin friction Ty at the interface is given by 


ae U,$4(0), 


In this case 


5092.13 
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The first approximation when py is large is 
C,, = $(1—A)(1+-A)}, (36) 
being in fact the limit as pu tends to infinity, which is the value for a solid 
lower boundary moving with velocity U,. 
3.2. A better approximation should be obtained by taking the functions 
¢, and ¢, to be quartic polynomials in yf and 7}. 
thus $y = + (1—e)(2nt—2n79+ 77) (37) 
and dy = c+ (A—c)(2nF— 238+ $4) (38) 
satisfy the conditions 
d(I)=1, (1I)=90, (1) =9, 
$,(1) =A, d,(1) = 0, $3(1) = 0, 
$;(0) = (0) = ¢, 


;(0) = $3(0) = 0, 


and 54 (0) - uhs(0), 
; 5+ pA 
provided that c = ~ P" as before (equation (32)). 
O-+p 
Then ' LF = 315 + 830€ eC” 
and I, = $A? t aoc = feec*, 


and the equation corresponding to (34) is 
6\? 7 ] [189(5 pe) 74+ 115A] 
ye py | (8/)(74A4-115)+ 189A] 


When pu is large, the first approximation to 8/p is 


(39) 


\ 741115A\1 
re) 74+-115A\3 f a>e 
pe 189A pu 
5 f a 
or —_= 0-857(up) t+ if A=O0. 
pe 


The velocity at the interface is given as before by equation (32), and 
the skin friction coefficient is given by equation (35), the first approxima- 
tion when py» is large being 

c 0-343(1—A)(1+-1-556A)!. 


To 
3.3. Good results have also been obtained in the case of flow along a 
flat plate by taking the velocity profile to be part of a sine curve. This 


suggests trying , : * 
o© ~— d, = c+(1—c)sin $anjz (40) 


and d, = c+(A—c)sin $75, (41) 
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which satisfy the conditions 
$,(1) = 1, $,(1) = 0, 
¢,(1) A, (1) 0. 
$,(0) = ¢,(0) = c, 
$;(0) = ¢,(0) = 0 

and 56,(0) ud,(0), 


provided that ¢ is given as before by equation (32). 


The 9 , , 
no Etenft3)roft-9) 


and the equation for 6/ is 
(°) l [ (2a 4)(8/u)+-(32—8)A+ (4 77) | 
U py [(8 u){(327—8)+(4—7)A}+4+ (27 4)A] 


The first approximation to 6/p Is 


( 1-425A+0-858\1 . 
; . \ if ADO, 
iv 2-283Appy 
) _ 
or 0-844(up)-* if A=O. 


b 

The velocity at the interface and the skin friction coefficient are given 
as before by equations (32) and (35), the first approximation to the latter 
when py is large being 

c,, = 0°328(1—A)(1+-1-657A)}. 

The velocity and skin friction coefficients at the interface have been 
calculated, using equations (39) and (42), for the values 5-965 x 104, 100, 
10, and 1 for pu, with A 0, and also for pu = 1, with A = 0-25, 0-501, 
and 0:75. They are compared in Table I of section 5 with the values 
obtained by accurate numerical integration of the differential equation. 
The agreement is quite good in all cases; the method of section 3.3 gives 
the best results for the velocity at the interface, and the method of section 


3.2 for the skin friction. 


4. The integration of the exact boundary layer equations 

4.1. Asymptotic expansions 

Before it is possible to start the numerical integration of the boundary- 
layer equations it is first necessary to investigate the asymptotic forms 
which the solutions of the equation 


2f"(n)+F(nF"(n) = 0 (43) 


take when 7 is large. 
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This equation has the property (see e.g. 1) that, if f(7) is any solution, 
then g(€) is also a solution, provided that 
€=an+b 

and f(n) = ag(€), 
where a and 0b are any constants. Then 

f'(n) = a9'(8), 

f"(n) = a9"), 
and so on, where dashes denote derivatives. 


We will first suppose that f’(1) > A as 7 > —oo. There are two cases to 
be considered, A = 0 and A > 0. In the first case, when A = 0, we can 
suppose that f(») > —a as n > —oo, where aisa constant. If we substitute 

a for f in equation (43) and integrate three times, we get 


f(n) ~ —a-+Aeim, 
where A is a constant. This suggests that we try the expansion 
f(n) A,+A,e#7+ A,e"-+ Agei7+.... 


Substituting in (43), and equating to zero the coefficients of successive 
powers of e#“”, we obtain the recurrence relations 


aA,+A,A, = 0, 
8aA,+4A,A,+ A}? = 0, 


27aA,+9A,4A,+5A, A, = 0, 


Hence A, a and A, may be chosen arbitrarily. A standard solution 
g(€) may be obtained by putting a 1, A, 1; then 
¢ . i. wee 7 oft | 
g(€) 1+et§— lef Seif izes (44) 


Any other solution with the same boundary condition can be obtained 
from this by making use of the property mentioned above. The expansion 
(44) is in fact convergent for € < 0. 

In the second case, if A > 0, we can assume that 


lim (f—An) B, 


n 
where B& is a constant. Equation (43) then becomes approximately 
1 PI : 

2f" +(An+ B)f” = 0, 


the solution of which may be written 


i B\2 
ae Anexp| (1 + 1) \ 
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where A is a constant, so that 
“5 lls 3 
f’ ~A4+ 2A erf | — I (y ! x) 


where erf z e dt, 


7 


und f ~ An+ B+ 2AM erfy n'(1 a 1} dy. 


Now it is known that an asymptotic expansion for erfz is 
l ¢ l 13 1 tas I 


eri z ~ l ; _ ° ° ° ~—-t 
) on | )~2 o> ® «ad » » » ~6 


Integration of this expansion leads to 


f~An+B a(1 ewe w); 


2>2 4-4 


where Z n(n } x): (45) 


A second approximation can be obtained by substituting the expression 
Ac ae" : ; 
\An+B __} for f in equation (43). We get eventually, when 7 is large 


und negativ e 


.. 3. O15 [Lae | _ 
par 2 (! 222 ' 474 -) SS  * — 
toe 1a . 
Jj ~ A 2AM erkz—- oA > a one (47) 
<. e . ¢ 22? 
and f” ~ Ade" — 4AM — (48) 


where z is given by (45). 
The corresponding asymptotic expansions for f(7) when 7 is large and 
positive, and when lim f’‘(7) 1, lim (f—7) = D, can be obtained in a 


n> n> « 

similar way, or from (46), (47), and (48) by using the fact that —f(— 7) 
is also a solution of (43), putting A = 1 and writing —C for A, —D for B. 
They are 


Ce- 3 15 1 C2e-2u* 
~ 7+D- l f a en ee 49) 
7 / yz Dap2 4u4 8 w ; ( ) 
; 1 C%e-2 
' ~w 1—2Cerfw 1 6 50 
J 4 w*' —" 
, 1a~ 

and f" ~ Ce ‘2 + (51) 

where Ww (7 +-D) 


and C is a further arbitrary constant. 
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4.2. Methods of numerical integration 
For the case when the lower fluid is at rest, i.e. when A = 0, the property 

mentioned at the beginning of section 4.1 makes it extremely simple to 
obtain a numerical solution for given values of the density ratio p and the 
viscosity ratio ». The standard solution g(€), defined by equation (44), is 
first integrated numerically from the differential equation, using (44) as 
a starting-point, until the zero &, of g(€) is reached. This function g(&) can 
be used in all cases to give the solution for the lower fluid, by choosing 
appropriate values for the constants a and b. The solution g,(€), for € > &), 
is then defined by the conditions 

91(€o) = 9, 

91(Eo) g'(&o); 

91(Eo) = (up)*g” (Ep). 
This function is integrated numerically, starting with a series solution to 
give a few initial values, until a constant value of g}(€) is reached. The 
constants a and 6b are then given by 

] 
a - aE PRT, b é., 
ViG(%)} 

and it is easy to see that the functions f,(7,), fo(y2) defined by 


fa(n2) = ag(€), 


film) = agi), 
where E=anetb (€ < &) 
and E=an,+5 (€ > &), 


satisfy all the conditions of the problem. It is evident that in this case 
there is a single infinity of solutions, depending only on the parameter pu. 

When A is not zero, there is a double infinity of solutions, depending on 
the parameters pu and A; and for given values of these parameters it is not 


easy to find the appropriate solution. If, however, py is given, a set of 


solutions with different (initially unknown) values of A could be obtained 
by starting at large positive values of 7, with the functions defined by the 
asymptotic expressions (49) to (51), using a set of values for the arbitrary 
constant C’, and integrating numerically to large negative values of 7, to 
find the corresponding values of A, at the same time satisfying the correct 


boundary conditions at the interface. Alternatively, the methods of 


sections 3.2 and 3.3 could be used to obtain approximate values of the 
velocity and velocity gradient on either side of the interface for given 
values of pu and A. These would define solutions of the differential equation 
for positive 7, and negative 7, which could be integrated numerically, 
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starting at the origin, until the velocity became constant in both directions. 


























‘ty By choosing the scaling factor a in the same way as before an accurate 
to solution would thus be obtained with a value of A nearly equal to the given 
he value. This seems in fact the best method and has been used in the only 
, is example worked out. 
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) 


starting at é 3, and using five places of decimals, of which four are 
retained in the table. The values for € - 3 were calculated from the 
series (44). 

Complete solutions have been obtained, in the case when the lower fluid 


is at rest (A = 0), for pu = 5-965 x 104, corresponding to air flowing over 
water at a temperature of 10° C., and for pu 100, 10, and 1. (The last 


solution has evidently been obtained before in the course of the work 
described in (1), but numerical details were not given.) For the case 
pu = 1, the solution has also been calculated when the lower fluid is in 
motion, with A = 0-5. The numerical results are given in Tables III-VII, 

nd the velocity distributions are shown graphically in Fig. 2. 

The principal features of these results are summarized in Table I below, 
which gives the values of the velocity ratio u)/U, and the skin friction 
oefficient c,, at the interface for the five cases. u,/U, is plotted against 

for A = 0, in Fig. 3, and (u,—U,)/(U,— U,) is plotted against A, for 


l, in Fig. 4. It is evident that the value of c,, depends only on u)/U;, 


whatever the value of A; c,, is therefore plotted against u,/U, in Fig. 5. 
> 


Values obtained by the approximate methods of 3.2 and 3.3 are also 


included in the table for comparison. 
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THE EFFECT OF THE NON-UNIFORMITY OF THE 
STREAM ON THE AERODYNAMIC CHARACTERISTICS 
OF A MOVING AEROFOIL 


By E. E. JONES (The University, Nottingham) 


[Received 5 September 1950] 
SUMMARY 


The problem of the motion of a cylinder of aerofoil section in an incompressible 
inviscid fluid, the velocity of which is non-uniform, involves the application of the 
two-dimensional complex potential function theory to a class of boundary-value 
problem which has not yet been considered to any great extent. A functional solution 
of the potential equation is determined, regular in the space external to a simple 
closed curve, which takes certain values on the boundary of the curve and certain 
values at infinite distance from this curve. A general solution in the form of a 
complex potential function is found in terms of these given values and the coefficients 
determining the profile of the cylinder section. 

The forces and couple on the cylinder, which moves in a general manner, are 
determined by use of the Blasius formulae for a moving cylinder. The effect of a 
general type of vortex wake extending behind the cylinder is also discussed. 

The results are applied to determine the effect of the non-uniformity of the stream 
on the aerodynamic charaeteristics of a moving aerofoil having small thickness and 


camber. 


1. Introduction 


THE air streams usually encountered in actual flight are such that the 


stream velocity varies from point to point in space. Another example of 


such a non-uniform stream is that existing in a wind-tunnel, due to the 
pressure gradient along its axis producing the effect of a convergent 
stream. In practice this necessitates a correction to the lift and drag 
acting on the aerofoil placed in the tunnel, and methods of evaluating 
this have been given by Taylor (1) and Goldstein (2). Taylor, making 
use of general hydrodynamical concepts but without reference to any 
particular type of stream, deduced expressions for the lift and drag on a 
stationary body in terms of certain pressure derivatives. Goldstein, by 
use of a two-dimensional transformation, deduced the forces and couple 
acting on a cylinder of aerofoil section at rest in a stream possessing 
curvature, the results of which reduce to those deduced by Taylor in a 
particular case. 

The complex potential theory developed in this paper has been success- 
fully applied by Morris (3) to the case of a cylinder moving in a general 
manner in a stream at rest at infinity, and in a uniform stream, and 
expressions for the lift, drag, and couple acting on the cylinder were 


[Quart. Journ. Mech. and Applied Math., Vol. IV, Pt. 1 (1951)] 
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deduced. The non-uniformity of the stream is also likely to have con- 
siderable effect on these characteristics, and in this paper general expres- 
sions are deduced as the cylinder moves relative to the non-uniform stream. 


2. General results 

The cylinder moves in an inviscid incompressible fluid with velocity 
omponents uw, v along axes Ox, Oy fixed in a right section of the cylinder, 
nd has angular velocity w about an axis perpendicular to this section. 
If z (= a+iy) defines the position of any point on the contour of the 
ylinder section referred to these axes, then the component forces X, Y 
re given by the Blasius formulae 

. . acy\2 
Yiix ko | (w+iwz)(@—iwz) dZ—tp | (@- ie —Z) ae 


Cc 


nd the couple lis the real part of 
q ; — a i. ae oe 
7) (w—+-tw2z)(W—tw2Z )z dz Sp (i twee —- ] 2ds— 


. 
( ¢ 





The integrals are taken in a positive anti-clockwise sense round c, the 
contour of the cylinder section in the z-plane. As usual w = u+iv, and 
Qis the complex potential function for the fluid motion round the cylinder, 
its gradient defining the complex velocity of the fluid. 

The cylinder profile is defined as the curve 7 = 0 in the conformal 
transformation 


z= > a, e*-™, C = €+in, (3) 


which transforms the space outside the cylinder in the z-plane to the 


inside of a rectangle in the ¢-plane, where 7 = +00 is chosen to corre- 
spond to z 1 a0, 

If the cylinder has a sharp trailing edge, the above transformation 
ceases to be conformal at this point, and in general this leads to an infinite 
fluid velocity there. This difficulty is avoided by assuming the Joukowski 
condition for finite fluid velocity at the trailing edge, which determines the 
value of the circulation. However, when a trailing vortex wake forms 
behind the cylinder the condition that the motion is irrotational does not 
ipply in the region of the wake, and thus Bernoulli’s equation for the 


pressure cannot be used at points on the cylinder in the neighbourhood 
F 


5092.13 
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of the trailing edge. To avoid this difficulty it must be assumed that the 
trailing wake is confined to a narrow band from the trailing edge, and the 
contour integrals can then be taken from the trailing edge, round the 
contour, and back to the trailing edge. Thus in transforming the contour 
integrals in the z-plane to integrals in the ¢-plane, the limits for ¢ are 
from ¢, to ¢,—27, where {, is the value of ¢ at the trailing edge, and @, and 
¢ are real on the contour of the cylinder section. 

3. The general method of determining the complex potential function 
for a cylinder moving in a given stream has been discussed by Morris (4). 
The function Q can be divided into two parts Q, and Q,, which are defined 
in the following manner: Q, is the complex potential for the motion of the 
cylinder in a fluid at rest at infinity, and Q, is the complex potential for 
the cylinder at rest in a given stream. These two motions can be super- 
posed by the addition of Q, and Q, to give the required motion relative 
to the cylinder in the stream. 

If a circulation « exists round the cylinder, then Morris (4) has shown 


that (21) . 
L 5 . = a 
; ¢ ik. 1 —e-tS—Ss) 
Q, > Qi ents - ‘log! aoe - | ds, (4) 
n=l 2a ) 2a \ ] —et(S—€s) 
s(%) 
where , e ; 
Q) = a,t@—d,w—iwb,, 
, eae 
i. An41W iwb,, in > 2), 
x 
and : = 
and Oy 2 Gia. 
r=0 


In this result the effect of a trailing vortex wake, extending behind the 
cylinder, has been included. The trail, forming a surface of discontinuity, 
extends from 2, to z, in the z-plane, where z, is the position of the trailing 
edge, and k, is the strength per unit length of the vortex trail at a distance 
s along the trail from the trailing edge. 

To determine 2, we proceed in the following way. If Q, is the complex 
potential of the undisturbed stream in the absence of the cylinder, then 
Q, must be such that it reduces to Q, at infinity and its imaginary part 


at the boundary of the cylinder must be zero. Divide Q, into two parts, 
F,(¢) and F,(¢), where F,(¢) tends to zero at infinity, and F,(¢) diverges 
there; F,(¢) diverges at infinity, and F,(¢) tends to zero there. To forma 
function which is to be purely real on the boundary » = 0, we may either 
add the conjugate of the terms in Q, treating ¢ as real, or subtract these 
terms altogether. The first part F,({) has a conjugate F,(Z) which by 
definition diverges at infinity, while F,(¢) tends to zero there. Conse- 
quently the required function is 

Q, 
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The forces and couple on the cylinder, moving in a general manner in 
a stream at rest at infinity, with a vortex wake extending behind the 
cylinder, have been determined by Morris (5), and these results will be 
quoted wherever they occur. These were deduced by the application of 
the Blasius formulae to the complex potential Q,, and the following 
nvestigation deals with the change in the characteristics of the cylinder, 


lue to the inclusion of the term Q, in the complete complex potential. 
4, The motion of a cylinder in the stream defined by the un- 
4 
disturbed complex potential: Q, = > A, z”. 


Although the method is applicable to the stream defined by Q, as a 
polynomial in z of any finite degree, the rest of this investigation will be 
be referred to axes 


mfined to a poly nomial of the fourth degree. Let 0, 


0,x,, O, y, fixed in space, hereafter called the stream axes, so that 


() 
ts | 
where the A,,’s are in general complex constants. This stream has a 


complex velocity A, at the origin of coordinates, and the terms con- 


l 
taining A,, Az, A, contribute to form the curvature of the stream. Since 
the velocity at a great distance is large, the results, if they are to have 
ny practical application, can only be expected to give the approximate 
effect of local curvature of the stream at the cylinder. 

The potential Q 


can be referred to the axes Ox, Oy fixed in the cylinder 
section, hereafter called the cylinder axes. If these axes are inclined at 
ingle @ to the stream axes, and the origin O is at a position z, referred to 


the fixed origin O,, then , » 1 optf 
“J ~ T we . 


9 


where z is referred to the cylinder axes. The complex potential then takes 





the form — 
@) . 8.2m 
—— ont ot] 
n=1 
where _ : P 
% A, 2,+A,27+A323+ Ag, 
‘ “ s we ~3\,i0 
O71 (A, +2492,+3Ag 23+ 4Ag ape", 
Oo (A,+3A 2,4 6A 423) », (6) 
8, = (Az+4A42,)e8%, 
O4 A,e* 


The z-plane is mapped on to the ¢-plane by equation (3), hence by 


ey iluating 22. 23 and z*, and substituting in 22,, we have 








— 


68 y. E. JONES 


where 
iN 1 9S ! 
M, = 8, Ay +28y 4941+ 385 a) 4,(49+a,)+ 48, a)(a§a,+ 3a, a, a.+4), 
Ms = §,a3+ 36, a2 a, + 25, a5(2a)a,+-3a,), 
ms = §,a,+ 46,48 4,, 
m, = 5,a6, 
and f(e%¢) is a function of positive exponential functions. Hence in this case 
n=4 " 
nal Y ; 
Fit) = > m,e-™, 
n=1 


and thus from equation (5), 
n=4 yy nas ; 
Q, : > mM, € “nts | Zz Mm, enil, 
n=1 n=1 


The complete complex potential function defining the motion of the 
cylinder in the given non-uniform stream is thus given by 


8(21) 
oO ¥ a ¢ Tt a a 
K ik. 1—e-“5—S)) 
Q= 3% enit_ KS | me | adi, ds +- 
n=1 27 E 2a | 1—eté ay) 


n=4 n=4 
+ ¥m, e+ ¥ m, eri. (8) 
n=1 n=1 

The force on the cylinder 

In order to evaluate the force on the cylinder, it is sufficient to deter- 
mine the contribution due to Q,. The additional terms arising in the 
Blasius formula of equation (1) are 
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(iv) In order to evaluate the integral 


Se 
* (EQ\? * (0Q\? dz 
bp | (=) az 2p | —| =, 
Z Oz 7 C4 dz/d¢ 
P &—27 


we note that the space outside the cylinder and wake in the z-plane 
transforms on to the inside of the contour in the ¢-plane, given by 
(a) the real axis from Z 
(h 


(c) the line parallel to the imaginary axis € = &,, 


= - 
<7 tO &, 


the curve from ¢, to ¢,, corresponding to one side of the wake, 


(d) the line v 0, 

(e) the line parallel to the imaginary axis 
(f) the other side of the wake curve from (,—2z to (,—2z. 
The contour is indented at ¢, ¢,—2z, and at these points we choose 


eQ/ef = 0, so that the fluid velocity at z = z, is finite. Also the curve 
corresponding to the wake must be indented at every point on it, since 
poles exist all along the wake. 

The above integral taken round this complete contour is zero, since no 
poles oceur within the contour. Since the discontinuity due to the vortex 
wake does not extend beyond ¢, the integrals along the paths (c) and (e) 
of the contour cancel each other. It has also been shown by Morris (4) 
that the integrals along the paths (b) and (f) reduce to 
p | k(%,—iw%,+¢@) ds. 


It thus remains to calculate the integral along a line parallel to the €-axis 


at 7 +-oo. At 7 0, we have 
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in the form > d,, e-ns+ fie), 
where 
dy, = x*/4n?+-32m,(a; 6+ m,—iwb,)+ 18m,(a, 6+ m,—1b, w) + 


8m,(a, 0+ M,— ib, w)+- 2m, (a, G—A, w+ Mm, — 1), w), 





d, = 24m,(a,@+m,—ib, w)+ 12m,(a, 0+ m,— tb. w) + 
+4m,(a, Gd) w+-m,— 1b, w)+-iKm,/z, | 
d, = 16m,(a,@+m,—tb, w)+ 6m,(a, F—G,w+m,—ib, w)— (9 
€ . | be 
— mi + 2ixm,/z, | 
d, = 8m,(a,@—G,w+m,—ib, w)—4m, m+ 3ixkm,/z, | 
d, —4m5— 6m, m+ 4ixm,/z, | 
7 9 ( 2 > | 
d, —8m,m,—12m, ms, d, = —9m5— 16m, m,, | 
} 
9 2 an 2 
d, 24m, m4, d, = —16mj, 
and f(e’s) is a function of positive exponential functions. Also at » = +« 
we write , 
(dz/df)-1 = 2 ¥ ce", 
= — ut 
n=1 
where 
. l . 0 . lq2 . 9 2 
Cj Qo; Cy » ° C3 As 9; C4 =A, ao; 
: 9 2) [73 Iq 
C5 = (day G43) /A, Cg = 4(ay 45+, 43)/A9, 
c, = (5a2a,+4a,a3+ 6a) a,a,+a3)/a5, (9a) 
> 9 ) 9 2 9 4 
Cn (6a5 a,-+- 8a) a,4,+ 12a, a,a,4+- 6a,a5)/a, 
Cy = (Tuga,+10a2a,a,+-9a2 a}+ 16a2.a,a,+9a3.a,+12a,a,a3+-a4)/ae. 
Since the constant terms in the integrand are the only contributions to 
the value of the integral at 7 +00, then these occur in the expressions 
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then the second integral reduces to 
7 riz, nil, | » 
Pp | k, P) +1 £ a > 3 Pnié ie | ds. 


In order to deduce the force acting on the cylinder, the above results 
are added to those already deduced by Morris (5) for the motion of the 
cylinder in a stream otherwise at rest. 

When a vortex trail exists behind the cylinder, following Helmholtz, 


we have 


and then the force is given by the expression 


r 1X 7 p(2. Bw Aw 1Dw) 7pw(2 bw —Aw t ?Dw)-+- 
; n=8 n=4 
imp > c,d,-4 2mpw\aygm,— > NG y 41 My,| 
n 1 n 1 
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2177p) Ay My > na, ,m iwi, p | k(etss e'Ss) ds— 
n=1 8(2) 
8(2,) 
p | k,(Z,—G,—Gy e+-dg e's) ds—ipay | k,(C,et+—l,e%) ds 
8(Z) 
ay 7 il,\ 
p | he) 2 P hs Pr+i' tt OS 
where , 
1 AyGgt > Ny 41 Gn41 
n=1 
B Ay A, 
( S nb, b 
/ vt 
1 
D by dy > nb, An+t: 


n=1 
The expression for the force now has four extra terms, all of which involve 
the effect of the non-unifor mity of the stream. 

It is also of interest to consider the value of the force when the vortex 
trail effect is not considered. It is then determined by the expression 

Yt+ixX imp(2. Bw Aw+iDwo)- 7pw(2Bu— Aw i Dw) + 
n=8 
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where the term in « has been modified from that deduced by Morris (5) to 
allow for the fact that the trailing edge corresponds to { = &,. 

Using the relations w = 6, and w = z,e- the functions 7h, can be 
determined as functions of 5,, and therefore of A,, from equations (6) 
and (7). 


6. The couple on the cylinder 
The couple on the cylinder is similarly deduced from the Blasius formula 





for the moving cylinder, and the extra terms due to Q, occurring in 
equation (2) are the real part of 
= ~aQ , {Od . C2 
p | D —* — iwz —*)z dz—hp | 34. —4)z dz+- 
Zz c EL ot ot 
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{:—27 
the real part of which is zero, since €Q,/éf = €©,/6f on the contour. 


The real part of the integral involving the acceleration terms is 
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It now remains to calculate the extra terms occurring in the integral 
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Considering this integral taken round the same contour as for the force, 
since no poles occur within the contour, the integral round the complete 
contour is zero. As before, the integrals along the lines parallel to the 
y-axis at €, and €,—2z cancel each other, and along the sides of the vortex 
wake curves the integral has the value 


p | k.(2,—iwz,+@)z, ds. 
At 7 0 we have 


dz\- , “ op 
i i} Sa, wit | ¥ c, enis! 
dv — } n j 


is q, enil 
=0 


7 
where g > a,€ 
The terms which contribute to the integral are the constant terms in the 


integrand, and these are given by 
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; n=8 ” (” =4 sc n=4 oe ) 
imp LdnGn—p | hs) 2 Ine ™™— D Ine”! ds, (13) 
n=0 “as n=1 n=0 
r=4 
where In = D> TM,Gr—n- 
r nm 


The real part of equation (11), and equations (12) and (13), when added 
to those deduced by Morris (5), give the couple on the cylinder moving 
in the given non-uniform stream with an accompanying vortex wake, as 
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In the expression for the couple extra terms occur, giving the effect of the | ere. 





non-uniformity of the stream. at infi 
The couple on the cylinder, when the effect of the vortex trail is not | are th 
considered, is eylind 
[ = imp(u*B—w B)— Lap(wD+@D)-+ lirp( Di— Div) —zpCw t values 
e ee 
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7. (i) These general results can be shown to agree with those deduced | 
by Goldstein (2) when the cylinder is at rest, with no vortex trail behind 
. : 9 
the cylinder. In this case w = w = z, = 0, then equation (10) reduces to tail 
a . n 8 
Y+iX = inp > c,d,, (16) 
n=1 
and the couple given by: equation (15) reduces to 
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‘ a osm Sa 7 y | 
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If ee ok *: v Wher 
f squares and products of A, and A, are neglected, and A, assumed real, 
these results reduce to those given by Goldstein, noting that the trans- 
. . . . . e < | 
formation used is the inverse of that defined by equation (3). 
(ii) When the cylinder is at incidence a to the real stream axis 0,2, 9 A 
and is at rest with its centre at the origin, and if there is no circulation iT 
; : 1 
round the cylinder, equation (16) reduces to f 
ola 
F 4 ~~ J J im Bia 
Y+ix dimpay A, Ay(dy€ - Ay e3!%), F accel 
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Here A and H are the energy coefficients for motion in a stream at rest 
t infinity, S is the cross-sectional area of the cylinder, and ép/éx and ép/éy 
are the pressure gradients along the stream axes at the origin when the 
cylinder is removed from the stream. These pressure gradients have the 


values 


pA,(A,+ As), ot ~ipA,(A,— Ag). 


cp 


C2 cy 


These results for the lift and drag agree with those deduced by Taylor (1). 


8 The condition that eQ/cZ 0 at ¢,, where 2'(%) 0, is equiva- 
lent to the Joukowski condition for finite fluid velocity at the trailing edge. 
However, it is known that this condition gives the correct pressure distri- 
bution over the cylinder only when the incidence to the main stream is 
small, and when this is so it is possible to determine « and « uniquely. 


For the stream under discussion the Joukowski condition gives 


K . ‘ _" mn 4 _ ie 
i> NG. ents is nm... e~"est-4 ” nm, ents | 
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Nei) ~¥ ¥ S 
ek en tee-$) elroy) 
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where also K ‘ ds 0. 


When the effect of the trail is excluded, then we have simply 


n=4 n=4 
is nQ) « i> nm, e-"S+i SY nm, ere, (19) 
Daw how } ya n hon n 


=i ? l n=1 n=1 


9, Aerofoils with small thickness and camber 

The forces and couple acting on a cylinder in the presence or absence 
of a vortex trail have been expressed as functions of the velocities and 
ecelerations of the cylinder, the coefficients determining the profile of the 
cylinder section, the coefficients determining the non-uniformity of the 
stream, and the position coordinates of the cylinder relative to the stream. 
Hence if the path of the cylinder is known as it moves relative to the 
stream, the change in the aerodynamic characteristics of the cylinder can 
be calculated at all points of its path. 

Although the effect of the trailing vortex wake has been included in 
these expressions, the integrals involved are extremely complicated. This 
is due to the fact that in a non-uniform stream the vortex wake does not 
extend in a straight line behind the cylinder but assumes a curved path. 
Thus in the further investigation of the effect of the non-uniformity of the 
stream on the characteristics of particular types of cylinders the effect of 
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the trailing vortex wake is not included, as further research is necessary 
before this problem can be pursued. 

In order to include the effect of the thickness and camber of the cylinder 
of aerofoil section, the following Joukowski transformation is used: 


1 


it? | <9 


<a. 


The magnitude of the thickness and camber of any given aerofoil is deter- 


mined by the value of z, in the transformation, and for a thin aerofoil with | 


small camber it is possible to neglect squares of z). The trailing edge occurs 
at ¢, where zp tae-% — —b, 


or 


= 


approximately b = a—2Xy. 


As shown by Morris (6), the transformation can be expanded in powers 


of z,e*S/a, and the values of the constants occurring in the equations for 


the force and couple deduced. 

If we consider the motion of the cylinder in the particular type of stream 
given by , , 
© . A, = —U, A, = —pU, A, = A, = 0, 
the undisturbed complex potential of the stream is given by 

fo To2 
Q, = —Uz,—pU2?, 
and further », U are assumed to be real, with » small. 

The cylinder is given a small anticlockwise angular velocity w about an 
axis perpendicular to the section, and has a velocity V (> U) along the 
positive x,-axis. The cylinder is at small incidence « (= @) to this axis, 
and at any instant has a position coordinate z, = z. = 2, say, referred to 
the origin O,. 

In this example, since « is small, the circulation round the cylinder can 
be determined uniquely from equation (18), and has a value 

os 2(U+V)(asin a+ y, cos «)—4a2uU sin 2a+-4apU 2 sin 

a ye SIN ~-+ Yo COS a) BU sin 2a+-4ayU x sin «. 
Substituting the known values of the variables and constants in equation 
(10), the lift Z and the drag D opposing the motion, are given by 
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4aV{3(a—2a)sin 20+ Yo(1-+-3 cos 2a)}+ 8a3w cos a+ 


mp 
+ 4a?w(V +4U)cos*a+4(U +V)?(asin a+ yy cos x) 
+ 4a?nU (2V —U)sin 20+ 16aU(U-+ V)ux sin a, 
D , ; Bd ; 
= 12aV {y) cos a+ (a—2ay)sin x}sin a+ 8a*w sin a+ 
™p 


+ 2a*%w(V +4U)sin 20+ 8a2uU (2V — U)sin*a, 
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ssary | and from equation (15), the couple is given by 
r 
7p 


ili ta*(a+-x,)V sin a—6a*y, V cos a—6atw+ 4a3wU cos a+ 
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2a(a—2x,))(U+-V)* sin 2a+ 8a3uU?(sin 3a—sin «)+- 
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8a°uU V (2 cos 2a—1)+ 8a2U(U+V)ua sin 2 
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When vp is positive, the cylinder moves as if towards a source with 


eter 
with 


-Curs 


velocity V, and angular velocity w at small incidence a. The effect of the 
curvature and divergence of the stream is to increase the lift, drag, and 
the couple on the cylinder as compared with that for motion in a uniform 


stream of velocity U. The lift and couple increase uniformly as the cylinder 


moves towards the region of divergence, the drag remaining constant to 


the same order of approximation. 
wers 


s for 


| When vp is negative, the cylinder moves as if away from a sink, the lift, 
lrag, and couple decrease below their values occurring in a uniform stream, 


und as the cylinder moves away from the region of convergence the lift and 


‘eam . 6 -s 
‘M | the couple increase uniformly, the drag remaining constant. 


The stream defined by 
Q, —Uz,— ipUz?, 


an be similarly investigated to deduce the effect of asymmetric local 


stream curvature on the characteristics of the cylinder. 


In conclusion, the author wishes to thank Dr. R. M. Morris for her 


tan | interest in this paper, which forms part of a thesis for the Ph.D. degree 


the of the University of Wales. 
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THE BEHAVIOUR OF FRAMED STRUCTURES 
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SUMMARY 
When a structure is subjected to several loads, each of which may vary indepen- 
dently of the other loads between certain limits, it becomes necessary to consider 
the possibility that failure may occur by the establishment of a cycle of plastic 
deformation, even though no possible combination of the loads could cause collapse, 
If, for example, a flexural member was bent repeatedly first in one sense and then 


in the other sense, so that yield occurred at each reversal of stress, fracture would 


be expected to occur after a comparatively short time. Alternatively, the permanent 


deformations might continue to increase each time a particular sequence of loads 
was repeated, so that large deformations would be developed after a few repetitions 
of the loads. In such cases the structure must be designed so that after a few applica- 
tions of the various possible peak load combinations, a state of residual stress is 
reached which enables all further variations of the loads to be supported elastically. 
When this happens the structure is said to have shaken down. 

In this paper the behaviour of framed structures under variable loads is discussed, 
and a theorem is established which enables framed structures to be examined with 
a view to determining whether shake-down will occur under a given set of loads 
which may alternate between prescribed limits. 


Introduction 

IN recent years methods of plastic design have been developed for struc- 
tures of ductile material, such as mild steel, which are based on the calcula- 
tion of the load at which a structure collapses owing to excessive plastic 
deformation. Such methods of design may lead to substantial economies 
of material, such as have been discussed by Baker (1). In most investiga- 
tions into the problem of plastic design it has been assumed that each load 
increases steadily from zero to its maximum value, and that the loads bear 
a constant ratio to one another. It often happens, however, that a struc- 
ture is subjected to several loads, each of which may vary a large number 
of times between prescribed maximum and minimum values independently 
of the values of the other loads at the same time. Such conditions of loading 
might arise, for example, in the steel frame of a factory building, which 
might be stressed by a wind load which varied independently of other loads 
on the frame, such as those due to overhead travelling cranes. The possi- 
bility then arises that although no possible load combination could cause 
collapse, a cycle of plastic deformations might be established in the struc- 
ture, no matter how often and in what order the loads were applied. For 


(Quart. Journ. Mech, and Applied Math., Vol. IV, Pt. 1 (1951)] 
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nstance, a member might be bent back and forth repeatedly, so that yield 
weurred alternately in tension and compression in its outer fibres. This 
would soon lead to fracture of the member. Alternatively, the permanent 
leformations might increase by definite amounts each time a cycle of loads 
vas repeated, leading to failure of the structure owing to the building up 
flarge deflexions. It is therefore necessary that in such cases the structure 
should be proportioned in such a way that plastic flow eventually ceases, 
state of residual stress being reached which enables all further variations 
fthe loads to be supported in a purely elastic manner. When this happens, 
e structure is said to have shaken down (2). 
The present paper deals with framed structures with rigid joints which 
rry load by virtue of the resistance of the members to flexure, and its 
biect is to establish the following theorem governing their behaviour 


inder variable loads: 


THEOREM. /[f any particular system of residual moments exists which would 
able all further variations of the applied loads between their prescribed limits 
to be supported ina purely elastic manner, then the structure will shake down, 
though the actual system of residual moments existing in the structure when 
thas shaken down will not necessarily be the particular system which has 


en found. 


The shake-down theorem was first stated in a general form by H. Bleich 
3). However, Bleich only proved this theorem for the particular case of 
structures possessing not more than two redundancies. A general proof for 
the case of trusses with any number of redundant members was given by 
Melan (4) and the proof given in this paper for the case of framed structures 
follows along somewhat similar lines. 

The analysis is based on assumptions which do not represent the be- 

viour of actual materials of construction very closely. However, this 
nables the presentation to be greatly simplified, and it is felt that the main 

tures of the problem are demonstrated. 


\ssumptions 
It will be assumed that the relation between bending moment and curva- 
ture is of the type shown in Fig. 1. Thus the behaviour is assumed to be 
istic unless the magnitude of the bending moment is J/,,, in which case 
the curvature can increase indefinitely while the bending moment remains 
nstant. In practice this assumption is not borne out, for elastic behaviour 
ses when the most highly stressed fibres reach the yield point, whereas 
curvature does not increase indefinitely under constant bending 


oment until the entire cross-section has yielded at a somewhat higher 
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bending moment. For instance, a typical mild steel beam of C-section, 
when loaded from a stress-free condition, may be expected to yield ata 
bending moment which is only about 0-8M,. Nevertheless, it is felt that 
the analysis which follows illustrates the essential features of the problem. 


M 





Mpp—— > > 








— 





Fic. 1 


Considering any point A on an unloading line, it is clear that the curva- 
ture may be regarded as the sum of the permanent curvature «,, and the 
elastic curvature x, = M/EI, where EJ is the flexural rigidity of the beam 
as ordinarily defined. Thus the behaviour can be summarized as follows: 

Either M = Elx,, M| < WM,, kK, = 0 
or M = 0, |\M| = M, K, # 9, K, = 0, 
where the dots denote increments of the bending moments and curvatures. 
This notation for increments will be retained throughout the paper. 

It is also assumed that the members are so proportioned that conditions 
of instability do not occur. The conditions under which this assumption 
would be justified have recently been discussed by Baler (1), who also 
showed that in most cases the effects of shear and direct thrust are negli- 
gible. A further assumption is that the structure is only subjected to 


concentrated loads, applied at a finite number of points. In this case the 
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bending-moment diagram consists of a series of straight-line segments 
between sections which are termed the joints of the structure. Joints 
therefore occur at any section where a load is applied, or where a member 
is supported or joins other members. It is then assumed that yield can 
only occur at the joints in a structure. This is clearly true in the special 
ease in which the cross-sections of the members are uniform between the 
joints, and is also justified in many other cases (5). 

If yield occurs at a joint, the curvature can increase indefinitely, and this 
implies that the slope of the member could undergo a discrete change at 
the joint. If this happens, a plastic hinge is said to have formed. In prac- 
tice, of course, such an abrupt change of slope cannot occur, since the 
strains must remain finite, but large changes of slope do take place over 
short lengths of members at joints, and the idealization involved in the 
concept of a plastic hinge has been justified by numerous experiments (1). 


The shake-down theorem 

Let -4, denote the bending moment at a joint i of a member, calculated 
on the assumption that the frame behaves elastically everywhere, and 
suppose that when all possible combinations of the external loads between 
their prescribed limits are considered, it is found that the maximum and 
minimum values of -@, are .@™®* and .@™",_ Further, let m; denote a set 
of residual moments at the joints, satisfying only the conditions of equili- 
brium when the external loads are zero. The actual bending moment 1/; 
at a joint is defined to be m “#,. Then the shake-down theorem states 
that if any particular set of residual moments m; can be found which 


satisfies at every joint the conditions 


A nax mM; : (M,); 


? ”) 
v/, mil mM; (M,,); } 


the frame will shake down, although the residual moments existing in the 
frame after shake-down will not necessarily be the m,. 

Before presenting the proof of this theorem it is necessary to state the 
principle of virtual work in the following simplified form suitable for the 
present discussion: . = 
| M’x’ ds4 » M's, 0. (2) 

In this equation M’ denotes any system of bending moments in equili- 
brium with zero external load, and x’ denotes any system of curvatures 
which is compatible with rotations y’, of the plastic hinges at the joints. 


The integration extends over all lengths of the members between joints, 


and the summation includes all the joints in the structure. 
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To establish the shake-down theorem, consider the positive definite 
quantity #, defined by 
: > (m—m)? 
E : ds. (3) 
2EI 
In this expression m represents the actual system of residual moments in 
the frame at any particular stage during the loading process, and m™ is a 
particular system of residual moments which has been found to satisfy the 
inequalities (1) at every joint. The integration extends over all the mem- 
bers of the structure. 
The change in E due to small increments in the applied load is 
-(m—m)m 
: ds, 
El 


Consider now the application of the principle of virtual work. The 


(4) 


system of bending moments will be chosen as (m—m), this system being 
in equilibrium with zero applied load. For the system of compatible 
displacements the changes of curvature and plastic hinge rotation asso- 
ciated with the increments » of the residual moments in the unloaded 
frame are chosen. Then from equation (2) 
*(m—m)m 
El 


Combining equations (4) and (5), it is seen that 


ds | > (m mb, 0. (5) 


. 


E > (m—myyp,,. (6) 


) 


Suppose now that at any section of a beam adjacent to a joint, 


m,—m, < 0. 
From (1), m; < (M,),;—-@?"*. Hence m;4--@™* < (M,);. 
Now m,+-@@ — M'**, where M/'** is the maximum possible bending 


moment which could arise at the cross-section 7 when the residual bending 


moment is m;. Thus M'e* . (M,,);- It follows that if yield is occurring 


at this cross-section, yp, < 9. 
By a similar argument, it can be shown that if m;+-m,; > 0, then 
us 0, so that — \ 7 rf 
p : (m;—m,)p,, > 0. (7) 


Comparing this result with equation (6) it is seen that L < 0. This 
shows that E decreases if yielding takes place, but remains constant while 
the structure behaves elastically. Since F is a positive definite quantity, 
it follows that it must eventually become zero or settle down at some 
positive value and thereafter remain constant, so that shake-down will 
occur. In the latter case the system of residual moments existing in the 


structure in the shaken down condition will differ from the ™;. 
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‘ite It should perhaps be noted that the residual moments m,, referred to in 
the foregoing proof as the actual system of residual moments at any par- 
. 5 ticular stage of the loading process, may not represent the residual moments 

that would arise in the structure upon removal of all the loads at this 


stage, for it is possible that some yielding might take place as the loads 


In 

_ were removed. At any joint 7, where the loads would produce a bending 

the ; moment -@; if the structure were entirely elastic, and the actual bending 

mn moment is V/,, m, is defined to be W;—-.%,. This would only correspond to 

the residual moment at the joint i which would arise if all the loads 

were removed provided that such a removal did not involve further plastic 

flow. 

t In conclusion, it should be mentioned that if, by some means of analysis 

| based on the theorem just proved, it is found that a structure will shake 

The } down, it is not necessary to carry out a separate investigation to ensure 

_— that there is no permissible combination of loads which will cause collapse. 
ible | \ structure that fulfils the conditions (1) cannot collapse, as has been ¢ 
mais shown by Greenberg (6). ' 
ded ‘ 
Discussion ’ 
} The assumptions made as to the relation between bending moment and : 

curvature, as indicated in Fig. 1, do not represent the behaviour of actual 

materials of construction very closely, owing to the range of bending 


moment which always exists between first yield and the development of 
fy] 
I 


plasticity over the entire cross-section of a beam. However, it is 


believed that the principal features of the behaviour of framed structures 


re described by the analysis based on those simplifying assumptions, and 
t should be remarked that a shake-down theorem has been established 
ing | which takes into account the range of bending moment which exists 
ie between first yield and full plasticity in a beam (7). 
‘te It has been shown that shake-down will occur if any set of residual 
} bending moments m; can be found which satisfies the inequalities (1) at 


, every joint in the frame, together with the conditions of equilibrium under 
the yj 1 
ero applied load. General methods are available for examining systems 
f linear inequalities of this type (8), and the application of these methods 
This to the shake-down problem has been investigated (9). 
vhile | In conclusion, it should perhaps be remarked that in many problems 
tity, where alternating loads have to be considered, the design may be governed 
some by the fatigue limit of the material if the peak loads may be repeated a 
will very large number of times. The criterion of shake-down is intended to 
1 the be applicable in cases where relatively few applications of the peak loads 






are expected to oecur. 
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SUMMARY 


The compatibility conditions for large strains are shown to be sufficient. <A 
m is proved which states that a given initial configuration and a given strain 
| determine the corresponding final configuration to within an arbitrary rigid 


ovement of the whole medium. 


1. Introduction 

Most of the early work on the theory of elasticity has been restricted for 
reasons of mathematical simplicity to the case of infinitesimally small 
strains:t and most of the recent work has been devoted to removing this 
restriction. 

A field of strain in a continuous medium can be described relative to a 
given coordinate system by six functions of position, known as strain 
components; not every set of six functions of position can be so used, 
however, because the strain components necessarily satisfy a certain set of 
six partial differential equations, known as compatibility conditions. These 
conditions restrict the way in which the values of the strain components 
can vary from one particle to another; the strain components at any one 
particle are independent and can be given values at will. The conditions 
are purely kinematic in origin, and their existence can be traced to the 
fact that the medium is assumed to be continuous and to deform in a 

mtinuous manner. When the strain is infinitesimally small and the 
oordinate system is rectangular Cartesian, the form of the conditions is 
well known; in fact one method of calculating the distribution of stress in 
in elastic solid is to solve the differential equations for the stress variables 
which are obtained by eliminating the strain variables between the stress— 
strain relations and the compatibility conditions. 


\n infinitesimal strain is generally defined as a movement for which the space deriva- 
of the displacement functions in Cartesian coordinates are so small that their squares 
products may be neglected ; it is worth noting that this requires not only that the ratio 
hange of separation to separation of every pair of neighbouring particles shall be small, 
that the rotational movement of every material element shall be small; the transla- 
movement not, of course, restricted. 


Quart. Journ. Mech. and Applied Math., Vol. IV, Pt. 1 (1951)] 
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Recently Weissenberg (1) and, independently, Green and Zerna. (2), 
have shown that the fact that the Riemann curvature tensor field vanishes 
in Euclidean space can be used to derive the compatibility conditions in 
a form valid for large strains and for any coordinate system. In their 
general form the equations are non-linear in the strain components, and 
are so complicated that their use in most practical problems is likely to 
present formidable mathematical difficulties. Nevertheless, it seems worth 
while giving some discussion of the conditions. 

In this paper we shall be concerned with two rather formal questions. 
We shall prove that the compatibility conditions as derived by Weissen- 
berg are not only necessary but sufficient, and we shall prove a related 
theorem which states that whenever an initial configuration is given, a 
knowledge of the values of the strain components throughout the medium 
is sufficient to determine the corresponding final configuration to within 


an arbitrary translation and rotation of the medium as a rigid whole. 


2. The description of strain 


We cortsider two configurations of a continuous medium, one of which 


we shall call the initial configuration and the other the final configuration. 


We suppose not only that the medium is continuously distributed in space 


but that either configuration can be obtained from the other by a con- 
tinuous motion, in the sense given below. (This restriction is physically 
important, for it means that any two particles which are infinitesimally 
near to one another in one configuration are infinitesimally near to one 
another in the other configuration, so that movements which, for example, 
involve diffusion of certain particles through the rest of the medium are 
excluded from the discussion.) 

It is perhaps simplest in the first instance to suppose that the medium 
is of infinite extent, covering the whole of space; it is possible to see, how- 
ever, that the arguments given below hold good for any region of medium 
of finite extent provided that this region can be regarded as a continuous 
part of a medium of infinite extent which moves continuously in going 
from either configuration to the other. 

Let us take any coordinate system (not necessarily Cartesian) fixed in 
space, and let us describe the position of particles of the medium in terms 
of coordinates referring to this system. Let (a!,a?,a°) and (X!,X?,X°) 
denote the coordinates of the positions occupied by a typical particle in 
the initial and final configurations respectively. Since one and only one 


particle occupies one and only one position in a configuration, the func- 


tional relations ; 
X¢ 


F(x) (a 1, 2,3) (say) (2.1) 
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which exist between the two sets of coordinates must be one-valued and 


nossess a one-valued inverse 


f'(X) ( 1, 2,3) (say). (2.2) 


We have here used x as an abbreviation for z!, 2”, 2°, and X as an abbrevia- 
tion for X!, X*, X°. Since the medium is continuously distributed in space, 
the variables X take on continuous ranges of values. 

By requiring either configuration to be obtainable from the other by a 
mtinuous movement we mean not only that the functions F*, f* shall be 
ontinuous but that all their partial derivatives up to at least the third 
rder shall exist and be continuous. 

Let us now suppose that the initial configuration is given, so that any 
particle can be recognized by the values of x, and let us regard the final 
nfiguration as capable of being varied. It is clear that a knowledge of 
he functions F’ in (2.1) suffices to determine the final configuration com- 


tely. In particular, the final configuration of any infinitesimally small 
sion of the medium is determined by the functions F’. The movement 

such a region from one configuration to the other can, however, be 
esolved into a translation and rotation involving no change in separation 
f particles of the region and a strain, which does involve such changes; 


e possibility of performing this resolution arises from the differentiability 


fthe functions #’. When stresses are set up in a moving material, they 
re primarily related to changes in separation of the particles and 
not directly affected by rigid movements, so it is convenient to have 
ie mathematical expression which describes the strain alone; it is clear 
it the functions F’, for example, will not serve this purpose. There is 
ict an unlimited number of possible expressions for describing strain, 


d several are in current use: 


t is generally the case that some are 
thematically more convenient than others in a given context; the form 
{the compatibility conditions will, of course, depend on the expression 
osen. We shall choose an expression which bears a close resemblance 
the compone nts of a metrical isor field and which in consequence 
ows us to derive the compatibiity conditions very readily by making 
se of certain known results: we shall indicate how the forms of the 
mpatibility conditions appropriate to other expressions can be obtained. 
Let g,,(a) denote the components of the metrical tensor field referred to 
the given coordinate system. The initial and final separations of a pair 


{ neighbouring 


particles are then given by the expressions 
q; (x) da'dx’, (2.3) 


) 


dS)? g;(X) dX*dX)!, (2.4) 
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where x, 2+dz are the initial positions, and X, X+dX are the final 
positions. Here and elsewhere we use the summation convention with 
dummy suffixes taking values 1, 2, and 3. 

We may now use (2.1) to express X in terms of x; substituting the result 


in (2.4) we obtain the final separation expressed in terms of the initial 


coordinates, viz. (dS)? = h,,(w) deride’, (2.5) 
oF eF8 

, TA . § Wiv>\\ Ia 

where hj (a ) J xB F(x)' at eal (2.6) 


It is clear from (2.5) that, the initial configuration being given, a know- 
ledge of the functions h;,(x) suffices to determine the separations of all 
pairs of neighbouring particles in the final configuration, and, moreover, 
that if the final configuration be varied, the values of h;,(a) change only 
when at least some of the separations change. Thus the functions h,,(x) 
may be used to describe the strain. (We note that h;;(x) are the components 
of a symmetric covariant tensor field, and so also are h,;(x) g;;(7); when 
h;; 


initial to the final configuration. ) 


(a) 9i;(*), the whole medium has moved rigidly in going from the 


3. The compatibility conditions 

In going over from the functions F'(x) to the functions h;;(2) we have 
achieved the desired result of finding expressions which are affected only 
by changes in separation of particles of the medium; we note, however, 
that the number of apparently independent functions has increased from 
three to six. This increase can only be apparent and not real because the 
six functions h,,(x) are completely determined once the three functions 
F‘(x) are known, the precise relations between them being given in equa- 
tions (2.6) above (the g,;(~) may be regarded as given once and for all with 
the coordinate system). By eliminating the functions F’ between these 
equations and such of their derived equations as may be necessary we 
should expect to be able to obtain certain equations for the h;,(x) and their 
derivatives which will effectively decrease the ‘functional degrees of in- 
dependence’ of the h;;(x) from six to three. These are in fact the com- 
patibility conditions for the strain components h,,(x). 

To carry out the required elimination directly would be so tedious as 
to be hardly practicable. Fortunately, exactly the same mathematical 
problem occurs in Riemannian geometry (in connexion with the conditions 
which a given set of functions g; (7) must satisfy in order to be components 
of a metrical tensor field in a three-dimensional Euclidean space) and the 
solution is well known; we shall be content to assume the truth of the 


required geometrical results and to show that they can be applied in our 


circumstances, not only to derive the compatibility conditions, but to show 








that 1 
the i 
but a 


into | 
Th 
for ou 
whos 
TH 
liffer 
stant 
from 
onste 
Eis 
quad 
dime 
reduc 
value 
the t 
trans 
We 
cond: 
that 
(2.5) 
the « 
Spaci 
quad 
syste 
tions 
quad 
Henc 
struc 
sions 
the f 
€ 


Cx 

































COMPATIBILITY CONDITIONS FOR LARGE STRAINS 89 


final that these are sufficient and to prove the related theorem mentioned in 
with the introduction. From this point of view, the present paper is nothing 


ut a translation of certain results from language familiar to geometers 


sult nto language familiar to rheologists. 
litia The theorem we require is attributed to Christoffel, but it seems better 
9 6 for our purpose to quote it in the form given subsequently by Eisenhart (3), 


hose proof is readily accessible. 


THEOREM (3.1). A necessary and sufficient condition that a quadratic 


ca ferential form q;(x) dx'dxi (i 1, 2, 3) be reducible to a form with con- 
¢ all dant coefficients is that the components of the Riemann tensor field constructed 
a from the coefficients g;(x) vanish; the transformation involves six arbitrary 
nstants. 
O! 
h Eisenhart proves the theorem for the general case of an n-dimensional 
ware iadratic form; we have put n 3 corresponding to space being three- 
when limensional. An inspection of the text shows that ‘reducible’ means 
, the | reducible by a transformation of variables having the properties of one- 
luedness, differentiability, and reversibility which we have assigned to 
the transformation 2 > X F(x) in section 2 above; we shall call such 
; ransformations ‘allowable transformations’. 
lave , . . “Ley: 
: We shall now use the first part of the theorem to derive the compatibility 
only a . P . 
onditions for h; (a). Let us suppose that two configurations are given and 
ever, | . . . 
that h;,(~) are the corresponding strain components. Equations (2.4) and 
trom — . . r _: as . 
| 2.5) show that an allowable transformation « > X exists which transforms 
e the , , . . . , TT 
bee the quadratic form h;,(2) dx'dx’ into the quadratic form g,(X)dX‘dX?. 
‘tions . . ; y re3awi ‘ 
Space being Euclidean, we know that g,,(X)dX'dX? is reducible to a 
qua . ; a . : 
a juadratic form with constant coefficients (by a change of coordinate 
with ‘ - . . ° 
ystem, in fact). Since the result of performing two allowable transforma- 
these . , — . ‘ ‘ i i 
ns in succession is itself an allowable transformation, it follows that the 
‘'y we ° . ° ° . . a i 
| quadratic form h,,(a) dx'dz/ is reducible to a form with constant coefficients. 
their 1 . ° . 
oi Hence, by the theorem, the components of the Riemann tensor field con- 
of in Fl ade . . ys 
structed from the coefficients h,;,(x) necessarily vanish. Using the expres- 
con . ° . ° _— ° 
sions for these components in the form given by Eisenhart (4), we obtain 
following equations: 
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Here h°? are functions of the strain components h;;(«) defined by the 


equations ' Ea : 
P Ah = 8 (i; 75 = 1, 2, 3), (3.3) 


6; being the Kronecker symbol. 

The equations (3.2) are the compatibility conditions for the strain com- 
ponents /;(x) in the form given by Weissenberg. They are evidently 
non-linear relations between h; (x) and their first- and second-order partial 
derivatives. Only six of the equations are independent (reference 3, p. 21); 
an independent set is obtained by giving the indices (ju, v: j, /) the values 
(a--1, a2: 6+-1. 6--2) (|, 6 1, 2, 3), with the convention of regarding 
as equivalent two indices which differ by a multiple of three (5). 

[t is interesting to note that the metrical components g;;(7) do not occur 
in the equations, so they have the same form whether the coordinate 
system is Cartesian or not. 

The terms on the left-hand side of (3.2) are of the first degree in the 
second-order derivatives of h,;(x); the remaining terms are of the second 
degree in the first-order derivatives. When the coordinate system is rect- 
angular Cartesian and the strain is infinitesimally small, h; (x) = 6;;, and 
therefore the terms on the right-hand side are negligible compared with those 
on the left-hand side; the equations so obtained have the well-known form. 
The possibility of obtaining the compatibility conditions in these restricted 


circumstances in this way has been pointed out by E. Trefftz (6). 


That the equations (3.2) are sufficient as well as necessary also follows | 


from the first part of theorem (3.1); to prove this we shall show that any 
set of twice-differentiable functions h; (x) h(x) satisfying (3.2) can 
describe a state of strain or, more precisely, that such functions being 
given, two configurations exist which are related by an allowable trans- 
formation 2 —> X F(x), where the functions F‘(x) are related to h;,(z) 
by equations (2.6). 

Since the given functions h;;(2) satisfy (3.2), it follows from theorem (3.1) 
that h; (2) dv'dz/ is reducible to a quadratic form with constant coefficients. 
The metrical form g;,(X) dX‘dX/ is also reducible to a quadratic form with 
constant coefficients. Two quadratic forms with constant coefficients can 
clearly be transformed into each other by an allowable transformation. 
Combining the various transformations involved, it follows that an allow- 


able transformation a2 — X F(x) exists for which 
h ij(*) dx'dai g;;(X) dX id XJ, (3.4 


Since values can be assigned at will to either set of differentials in this 
equation, it follows on making the substitution X F(x) that the functions 


F'(x) and h,,(x) are related by equations (2.6). Further, since « > X = F() 
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san allowable transformation, we may take any continuous configuration 


fthe medium and regard it as the initial configuration and then use this 


transformation to define the final configuration. The given functions h; (x) 


then describe the corresponding strain. This completes the proof of the 


The form of the compatibility conditions given above is appropriate to 
use of the functions / s strain components; some other types of 


in components which have been used are expressible in terms of h; (7) 


the metrical components g;,(7), and the form of conditions for these 
n be obtained at once from the conditions given above; for example, if 
use h;,(2)—g;(v) = e;,(x) (say), we merely substitute e;(7)-+-g,(xv) for 
in (3.2). and thus obtain conditions for €;;(«). There are, however, 

in components which cannot be expressed in terms of h; (2) and g; (x) 

ne. These arise when the separations of particles are described in terms 
the coordinates of their final positions, instead of their initial positions 
sabove: the two methods of describing changes in separations which arise 
this way have been called (7) ‘Eulerian’ and ‘Lagrangian’ respectively. 
some purposes the distinction between Lagrangian and Eulerian 
scriptions is important, but this is not the case for the questions dis- 
ssed in the present paper; this follows from the fact that we consider 
two configurations which are on an equal footing; although we have 

en preference to one configuration by describing strain in terms of the 
dinates of initial positions, it can be seen that the whole treatment 
equally well be carried through in terms of the coordinates of final 
sitions. Thus if we introduce quantities H, (X) such that H, (X) dX'dX/ 
the square of the initial separation of particles finally at X, X+dX, a 
ple interchange of the words ‘initial’ and ‘final’ wherever they are 
ied to configurations in the arguments given above allows one to 
uce that the compatibility conditions for H;,(X) have exactly the 
e form as the conditions given for h;,( r) in (3.2) above. Other Eulerian 
pes of strain components are expressible in terms of /;,(X) and g;(X), 


the compatibility conditions are accordingly readily obtainable. 


+, A property of the strain field 
Let us suppose that a coordinate system and an initial configuration are 
en, so that the coordinates of the initial position of every particle of the 
edium are known. We shall now consider the extent to which a know- 
ge of the strain components /;,(x) allows us to determine the correspond- 
g final configuration. Because this information allows us to determine 
final separation of every pair of neighbouring particles, we should 


pect the final configuration to be determined to within a movement 
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which involves no change in the separation of neighbouring particles 
and therefore no change in the separation of any particles. The second 
part of the theorem quoted above enables us to show that this is in facet 
the case. 

We have the same conditions as in the proof of the sufficiency of the 


compatibility conditions set out above where, given quantities h(x) } 


satisfying (3.2) and an initial configuration, we established the existence of 
a final configuration for which h;,,(x) dx‘dx’ gave the square of the final 
separation of particles initially at (x), (w+dx); we now consider the 
uniqueness of this final configuration. Let us first consider the case in 
which the given coordinate system is rectangular Cartesian. Let X,X-+dX 
be the coordinates of the final positions of particles initially at 2, 2+dy, 
We then have h,,(x) dx'dx’ = > (dX")?, where x + X is an allowable trans- 
formation; the second part of Theorem (3.1) applies, and shows that this 
transformation can involve at most six arbitrary constants. Consider the 
variables X‘, defined in terms of the variables X‘ by the equations 


Xt = At+ BiXi (i = 1, 2, 3), (4.1) 


where A‘ and B} are constants, and B} are the components of an orthogonal 
matrix. Clearly 2 —> X is an allowable transformation, and 


h,,(x) dxtdzi = } (dX')? = > (dX*)?, 


so that 2 > X defines another final configuration consistent with the given 
data. It is evident from (4.1) that this final configuration is obtainable 
from the other final configuration by a translation and rotation of the 
medium as a rigid whole. Moreover, since three of the constants B} can be 
given values at will, the transformation X — X and therefore x > X 
involves six arbitrary constants, so that all final configurations consistent 
with the given data are obtainable in this way. Hence the strain com- 
ponents h,,(~2) determine the final configuration to within an arbitrary 
movement of the medium as a rigid whole. 

It is easily seen that the same result holds good when the coordinate 
system is not rectangular Cartesian. From the given data, h;;(2) da'dz’ is 
a scalar for arbitrary values of the differentials, and therefore h,;(x) are 


components of a tensor field; the compatibility conditions express the 
} I : ) 


vanishing of the components of a tensor field; therefore the strain com- 
ponents induced in any other coordinate system by the given strain 
components h,,(x) will also satisfy compatibility conditions of exactly the 
same form as those satisfied by h,,(x). In particular, we may transform 
to a rectangular Cartesian coordinate system, and the arguments given in 
the preceding paragraph can then be applied. 
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CONTACT STRESS AND DEFORMATION IN A THIN 
ELASTIC LAYER 


By MARGARET HANNAH 
(Dept. of Textile Industries, University of Leeds) 


[Received 4 July 1950] 


SUMMARY 

The rollers used in drafting textile fibres are disk-shaped, the lower one made of steel, 
the upper with a metal centre and a fairly thin elastic cover. The length over which 
the rollers are in contact, and the pressure distribution over this length, are factors 
which affect their performance as drafting agents. The purpose of this paper is to 
show how these quantities vary with the material and thickness of the cover, the 
pressure between the rollers, and the roller size. The effect of allowing slippage at 
the inner boundary is also considered. 

The system can be considered mathematically as one of generalized plane stress 
in an elastie layer, with given displacement conditions on its inner boundary (the 
interface between metal and cover will be termed the ‘inner boundary’ of the cover 
and subject to pressure by a body of given shape on its free face. The layer is 
sufficiently thin for the ifner boundary conditions to affect the stresses in thi 
contact zone. 

The analysis of contact stresses was first carried out by Hertz (1), and quoted in 
Love (2). The application to the two-dimensional case is given by Thomas and 
Hoersch (3)—their results, which are for plane strain, may be converted by the usual 
modification of the Poisson’s ratio to those of generalized plane stress. These analyses, 
however, only hold if the contact stresses are the only effective forces over the contact 
ZOTLE 

The effect of the boundary conditions on the solution for a single isolated force 
may be found straightforwardly by a method given in Coker and Filon (4). The 
displacement due to any pressure distribution over the contact zone can then be 
determined, and the actual pressure distribution may be found by imposing thi 
condition of known displacement over this zone. It has been found most convenient, 
in practice, to determine the difference between the pressure distributions for an 
infinite and finite thickness ; this difference can be expressed as a Fourier series, and 


a sufficient number of the coefficients can be found to give any desired accuracy. 


Solution for infinite thickness 

Ir the free face is the plane y = 0 and the y-axis is measured vertically 
downwards, then it is known that the displacement v in a semi-infinite 
plate due to a downward force W at the origin is given by 


9 
- ———_ logr+ constant, (1) 


9 


where £ is the Young’s modulus of the material, and 7? = x?+-y’. 


[Quart. Journ. Mech. and Applied Math., Vol. IV, Pt. 1 (1951)] 
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It can be shown from this that, when an elastic disk of diameter D, is 
pressed against a hard disk of diameter D, with a force W, contact is over 
length 2h, where 
2WD 
h? i and 1/D 1/D,+1 D,. (2) 
tis assumed in the derivation that & is small in comparison with D, and D,. 


The pressure distribution over this contact length is given by 


21] 
P r) 2| | v2 hh?) ( h M h), (3) 


ath 
where P(x) is the force per unit length at x. These results may be derived 
from those given by Thomas and Hoersch (3) (with a slightly different 
notation 

The displacement function for the idealized (and non-physical) isolated 
force contains an inevitable logarithmic term giving a non-physical infinite 
lisplacement at the origin. However, the physical separation G(x) between 
the two surfaces at points outside but near to the contact zone may be 
found. We have 


G(x) i x? / D—v(h)—h?/D, 
e GQ’ (a) v’(a 22/D and Gh) 0. (+) 


By using (3) and (1) and integrating over the contact zone, 


t) - (1 xv” h?) log(a—a ‘) dx’, 
a2hh da 
which, from (2 ecomes 
2h/ Dfx/h—(a?/h?—1)8} for |z h, 
nd so (7(2) 2h D | x2 he 1) da 


-cosh-lx/h'. 


























96 





M. HANNAH 


Isolated force on a thin elastic layer 
A solution for the stress function following the methods of Coker and 
Filon (4) with some slight changes of notation, is 


rf 


* 
| COS 2 


X = —|[{o1(m)-4 


. my $o(m)}cosh my+-{44(m)+ my d,(m)}sinh m y| dm, 
2m ‘ ; ia ; ‘ 


0 (6) 
which can be made to satisfy four sets of boundary conditions. 
Coker and Filon give the solution appropriate to an isolated force W 
on the free face, with an inner boundary in contact with a smooth plane 
so that there is no shear and no vertical displacement there. This solution 


f 2W , 4W\  sinhw 
Pim) 7 Pam) T = sinh 2u’ 
4W sinh2x 4W cosh u 
d.(m) = ; d,(n : : : (7 
7 ( T is +-sinh 2u aim) ( 7 lia +-sinh 2u | 


where u mb, and 6 is the thickness of the layer. A constant term, 


2W cr 1 
2 dm 
» 2 
7 m am* 
0 


must be subtracted to give convergence at m = 0. The corresponding 
vertical displacement of points on the free boundary is 


ra 


vials 2W ff (coshz—1)ecos zax/2b d (8) 
sta ; * . 
vTE . 2(z+-sinh z) 
0 


where z = 2uw. 

Proceeding on similar lines, the solution for an isolated force W and an 
inner boundary having no horizontal or vertical displacement has been 
found by the present writer to be 


$,(m) = 2W/z, 


$o(m) = —$g(m) = | 


2W\ 
4(m) == ( che L 


7/(1+ 


2W\ {(3—4e)cosh usinh w—u?' 
u 


2+. (1—2o)?+ (3—4c)cosh?w’ 


7 


(1— 204 


(1—2o0)?-+ 


cosh?x) (9) 


(3—4c)cosh2u’ 


where o n), n being Poisson’s ratio. The same constant term as in 
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2W\ ¢ 1 ; 
the above solution | —, dm, must be subtracted to give con- 
Tv -m- 


vergence. The corresponding vertical displacement of points on the free 
boundary is 


2W fy (3—4o)sinhz—z [= dz 
. : , COs : 
ah ; | 12° ) 4c0)cosh z+ 5— 120+ 8c? | 2b, Z 


(10) 


These solutions are for a straight layer on a horizontal plane. They may, 
owever, be applied to an annular layer if the thickness is small relative 
to the diameter, i.e. if b is small relative to D,. 


Determination of contact stresses 
Equations (8) and (10) give the displacements for an isolated force with 


two different sets of boundary conditions, in the form 


ow ze 
(a) i | Q(2)e0s( 55) dz. 


0 
The displacement due to a pressure distribution P(x) over —h <a <his 


x 


2 Of pn ae [ ar. f**—2 Ni 2. 
r) 7 P(x’) da | Q(2)cos) h io 


« 


0 
nd it is required to determine P(a) so that this shall represent the given 
lisplacement, d—a?/D, over the contact zone —h <a <h, d being the 
deflexion at the midpoint of the zone. Because of the logarithmic infinity, 
the condition has to be differentiated with respect to x; such differentiation 


s permissible under the integral sign. Thus it becomes 


eo ._ (z(a—2’)) — 
27/D P(x") da 2()(z)sin dz h<2z<h), 
ih | an ( ) 
or, inverting the order ot integration. 
9x/D Q(z) dz , P\ inl ©) ay’ (—h<2<h). (11) 
—s Cd | 2b | 
\ssume now, by analogy with (3), that 
“ 2n—1)an 
P Yr) wih (] v2 /h?) -— S a, COS 1 pa (12) 
| ho 2h J 


n=1 
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an even function, which makes the pressure zero at x th. Then 


h 


r= [ P(r) dx = whl jet SY (—11{_44_| 
W Hg P(x) da wh | bo > | - \(2n—1)xJ | 


h w= 
And 
a 
| Pta'ysin| _ } da 
h 


1)cos(zh/2b) 
(hz/2b)?—{(2n—1)37}? 


La 


n=1 


2me| (I 2)oF,(zh/2b)+-(h/2b) YS (—1)n un 


where J,(x) is the Bessel function of the first order. Writing h/2b 
xh = X, (2n—1)}7 = yw, the relationship (11) becomes 


va 


X /D = (w/Eh) [ 2Q@)[ 0x2) z+K >( 1) 


, 4,(2n—1)cos Kz) 
y Kz? y?? | 
0 

for —1< X < I. 


Putting 


(—1)" (2n - ] pee v F (zx) 


x of, 


(14) becomes 


x 


X/D = (w/Eh) { 2Q(z)| F(Kz)+ 3, F,(Kz)|KsinKXzdz (—1<X; 


0 n=1 


And we have also, from equation (4), 





, | sin za/2b, 


K, 


sin K Xz dz 


(14) 


(15) 


A}. 


(16) 


(1/h)G’(a) = —(2w/ Eh) ( 2Q(z)| F(Kz)+ > a, F,(Kz)| K sin K Xz dz+-2X/D 
. 1 


0 ss 


(|X| > 1). 


(17) 


By inspection of (8) and (10) zQ(z) > 1 as z > o0 in both cases, so that it 


is convenient to consider first 


x 


H,(X) = { F(K2)K sin KXzdz 


0 


; 


and H(A) = [ F (Kz)K sin KXz i.| 


0 


(18) 
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These are known integrals, whose values are 


H,(X) X (—1<X <1) 
and X —(X?—1)! (|X| > 1) (19) 
W(X) 1)"-! { cos(l1—X)xs,, Si(] X )ys,, —cos(1-4 X xb, Si(1 +X )u,, ) 
7 \—sin(1—X)b,, Ci(1—X)h,,+sin(1+ X)h, Ci(1+X)ph,,]’ 
(20) 
where Si(x) wasted | Ci(2) = du, 
J u d u 
is usual. 


From (19) it can be seen that, as b, the layer thickness, tends to infinity, 
condition (14) becomes 


X/D = (w/Eh)\X+ Ya,H,(X)} (—1<X <1). 
1 


nt 
The solution of this and equation (13) is a, 0, all n, w = 2W/ah, and 
i? D = 2W/7E, which agrees with the previously obtained solution, as 


quoted in (2), (3), for this case. Define 


L,(K, X) = | {1—2Q(z).F(Kz)K sin KXz dz, 
L,(K, X) = | {1—2Q(z)F,(Kz)K sin KXz dz. 


| (21) 
| 


Then (16) becomes 


X/D = (w/Eh)|H,(X)—L,(K, X)+ > a,{H,(X) L,(K, X)}] 
n=1 


(_1<X <1). (22) 

Equation (22) can now be used with (13) to determine values of w, h, 

and a, for any given K. The a, are first determined by eliminating 
wD/ Eh); then using these 


, H,(1)—L,(K,1)+ > a,{,01)—L,(K,1)}] aw 
h —" eed 1 (23) 
D lat ¥ (—)"“{4a, /(2n—1)z} | = 
l 


And when K 0, from (2), 
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Thus, for the same loading, the contact length is smaller for a thin layer 


in the ratio 


CK) hi/ho, (24) 
Hy(1)—L)(K.1)+ > a,{H,(1)—L,(K, 1)}]3 
here C,(K — , 25 
where (A) " y “yf, |(2n—1)x2) (25) 
n 1 


Similarly the loading required to give the same contact length is given by 

Wi. = CyT*M. (26) 
The contact pressure at any point, P(X), may be compared with the 
average pressure over the zone, W'/2h, to give a non-dimensional function 
R(X); when K 0 this function is (4/77)(1—X?)!. From (12) and (13), its 
value in the general case is 


2h P(X) 21(1— X*)*-+4 ¥ a,, cos }(2n 1)Xz| 


R(X) = : oan 
HW lat ¥ (—)"“ 4a, /(2n—1)z} 
n=1 
so that R(X) = C,(K, X)R,(X), (27) 
(1—X?)!+ a, cos }(2n—1)Xa 
where CK, X) > w=) ee (28) 
1+ > (—)"“{8a,,/(2n— 1)a?}}(1— X?)! 
n=1 


C,(K,0) will be written simply as C,(4) in what follows. The ratio of 


the actual contact pressures at X 0, the midpoint of the zone, for the 
same contact length 2h, is then 
(3 CU, Ci. (29) 
Further, from (17), using (18), (19), and (21), we find 
G(X) 
2uht. a ; : or) a ‘ — 2Xh? 
z [x-(x 1)! LK. X)+ ¥ a,{H,(X)—L, (KX )}] + 5 
(30) 
Write 
x x 
| H,(X)dX =1,(X) (n > 1); | (X?—1)§' da = N(X) (31) 
i x i 
| L,(K,X) dX M,(K,X), all n. 
i 
From (5), 
x 
£G(X)} 9 = (2h2/D) | (X2—1)! dX = (2h2/D)N,(X). (32) 
1 


Then, if the loading is such that the contact lengths are the same, using 


eh Sa G(X) = CK, X)G)(X), (33) 
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CONTACT STRESS AND DEFORMATION 
where C,(K, X) Re l Gs )(x9) (54) 
(2) = SN 1+ 3a, | 
n=1 ; 
wnd 


A(X) = X2—1—2N,(X)—2M,(K, X)+2 3 a,{1,(X)—M,(K, X)}. (35) 


nt 
n=1 


It is found that, using the term in a, only, equation (22) can be satisfied 


to within 3 per cent. over the whole range of X for K = 1, which is the 
worst case. As a check on the adequacy of this approximation in the rest 
of the work, a, and a, were also obtained for K | by working out H,, 
H,, L,, L, at X (0-25, 0:75 and solving the four equations given by (22 
ut X 0-25, 0-5, 0-75, 1. The values found were a, 0-413, a, = 0-0 

fl, 0-002, as compared with a, = 0-511 from the first approximation. 


It was found that the values for C, and C, (equations (25) and (28)) for the 
two approximations differed only in the third decimal place, while C, 
equation (29)) differed by one in the second decimal. The results in 
Tables | and II were worked out on the more approximate formula, and 
so are given to two decimal places only, and the values for C,; have been 
corrected. The results should then be reliable to this order. Figures are 
given for C,, C,, and C, for a fixed inner boundary, » = 0, and a range of 
values of A; for comparison values have also been found, for A = 1 only, 
when 7» } and also when there is sideways slip at the inner boundary. 
Values of R(X) have also been calculated for X between 0 and 1, for 
K=0Oand K 1, and are compared in Fig. 1. 

C,, C,, and C, all refer to properties within the contact zone, while C, is 
1 measure of the deformation outside, but near to this zone. Here it is 
found that the first approximation is inadequate, when K = 1: the com- 
parative figures for C,(1, 2) are, from the first approximation, 1-242, and 
from the second, 1-165. The labour of extending this for the other 
values of A was not thought worth while for present purposes, since C, is 
the least important of the ratios calculated. The second value, 1-17, for 
(,(1, 2), may still be in error by two or three in the second decimal, and 
should be regarded as showing only the order of the effect. 

The solution appropriate to any given layer thickness and given loading 
is found by using (2) to get, for the given conditions, the contact length 2h, 
for an infinitely thick layer. Then the value of h for a finite thickness b 
must satisfy the two relations (14) and (24), 


h 2b kK, 


CK ho. 
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so that K is defined by the relation 
C,(K) = (2b/h,)K. (36) 


The other quantities can then be found from the table. 


Details of computation 

For small values of KX, L,(K,X), L,(K,X) are best evaluated by 
Simpson’s Rule in the usual way. If AX is large, however, the rapid 
oscillations of sin A Xz would require a very small length interval for z, 
and in this case it is better to use a different integration formula which 
will now be derived. 

By fitting a parabola to q(x) through the points q(n7z), g{(n+4)z}, 
q{(n+-1)7}, it can be shown that 


(n+1)7r 
F ; 7?2—8 a r7—8 . 
q(x)sin x dx = (—)" 5— 4(n7r)+-(16/7?)q{(n+-4)a} + ——— q{(n+ 1)z} 
nt : 
(—)"(0-189q,, + 1-624q, ,,+0°189¢,,.,) (37) 


with an error of 0-0562q'", q'Y being, as usual, the value of the fourth 
differential at some point in the range [nz, (n+ 1)z]. This is more accurate 
than the corresponding Simpson formula, with strip length 47, whose 
error is 0:1063q'". 
Any number of such strips can clearly be combined, and, because of the 

alternating sign, the end-values of each disappear, so that 

Na 

| q(x)sin a da (—)"[0-189q,,+ 1-624{q,, "dase tt 


nT 


H(—)8tgy y+ (—)"10-189qy]. 38)" 


In the computations for Ly, L,,, the integrand is of the form q(z)sin pz, 
where p = KX, and q(z) > 0 as z—> oo. Taking y as the largest value of z 
for which q(z) is appreciable, 
Y PY 
q(z)sin pz dz - q(z/p)sin z dz, 
) 


0 0 
which can be evaluated by the formula (38) using py/z strips, and so with 

a maximum cumulative error of 
0-056 2yq'* sae she 

YY — 0-0] 79yq'* p 1 
7p 

* It has been pointed out to me that this formula is a special case of a general formula 
derived by Filon (5) for evaluating integrals of this type. He deals with strips of width 6, 
and obtains an expression in terms of sums of odd and even ordinates of the function 
q(x)sin x; he gives numerical values of the coefficients which occur for values of @ between 
0 and 45°. Formula (38) above corresponds to @ 47 in this notation. It may still, however, 
be of interest in itself, in that it involves only sums of ordinates of g(x), and so is easier to 
evaluate than the more general case which brings in q(#)sin 2. 
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Similarly, using Simpson’s Rule with a length interval h, and an error per 
strip of g*h®/90, the cumulative error is 0-0056yq'*h*. Thus Simpson’s Rule 


is better only so long as 


p*h* " 3°218, 
e. ph < 1-34, 


In the computations h 1, and Simpson’s Rule was used for values of 


KX up to 1-2, replaced by formula (38) for larger values. 


Results of computations 
TABLE | 


0, inner boundary fixed 


7 
/ 
} 6 8 I 
788 211 0°357 O'sII 
1 I 5 r ) 73 ‘O7 
I 1°03 1°05 1:08 
1°63 1°97 2°37 
TABLE II 
K | 
il inner boundary 
1-08 — 
5 2°43 
Considering first the set of results for 7 0. it can be seen that the 


] 


effect of thickness on the deformation is quite considerable. The length 
| of the contact zone is reduced to nearly two-thirds for the same loading, 

is K increases from 0 to 1; or if the same contact length is produced, the 

total loading has to be more than doubled. The deformation outside the 
ontact zone is also considerably changed, so that the separation between 
the two surfaces, under conditions giving the same contact length 2h, is 
10-20 per cent. greater for the thin layer, K 1, than for the infinite 
layer, at a distance h from the nip. 

The pressure at the midpoint of the contact zone is very greatly increased 
when a thinner layer is loaded to give the same contact length. This is 
mainly because of the increase in total load, but the values given for C, 
show that a slightly larger proportion of the load is effective towards the 
centre of the contact zone, as the layer becomes thinner. This effect can 
be seen more clearly in Fig. 1, which compares the pressure distribution 
Oand AK = 1. 


over the contact zone for K 
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Poisson’s ratio can be seen to have little effect; the maximum possible 
change, from 7 = 0 to n = 3, produces, at K = 1, a 4 per cent. decrease 
in contact length, or an 8 per cent. increase in load for a given contact 
length, with a corresponding 9 per cent. increase in pressure at the mid- 
point of the contact zone. 

The results for a layer which can slip at the inner boundary lie between 
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those for 7 = 0 and 7 = }, and are very near to the former set. This is 


to be expected, since the freeing of the inner boundary obviously has a 
much smaller effect on a substance with less tendency to spread sideways. 


Conclusions 

Next to the roller diameter, and elastic modulus, the layer thickness is 
the most important factor determining the relation between loading and 
deformation for this type of roller. A greater load is necessary to produce 
the same contact length, and even then the gap between the two rollers 
near the contact zone is much larger, for a thinner elastic layer; this last 
property results from a smaller vertical deflexion, i.e. a smaller approach 


only slightly affected by layer thickness. 
A change in Poisson’s ratio makes an appreciable difference to the 
loading necessary for a given contact length, but has very little effect 








of the roller centres. The pressure distribution over the contact zone is 
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otherwise; slipping at the inner surface gives results very similar to those 


with fixed inner surface and a zero value of Poisson’s ratio. 
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THE FIELD INDUCED BY AN OSCILLATING 
MAGNETIC DIPOLE OUTSIDE A SEMI-INFINITE 
CONDUCTOR 
By A. N. GORDON (Imperial College, London) 
[Received 14 March 1950] 


SUMMARY 
Formulae are given for the external field produced by an oscillating magnetic 
dipole located outside a uniform semi-infinite conductor. The normal component of 
the field induced by a circular alternating current filament at the surface of the 


conductor is also considered. 


1. Introduction 
It has been suggested that information about the variation with depth of 
the earth’s conductivity can be obtained by means of a large-scale local 
experiment over a suitable flat region of the earth’s surface. In this ex- 
periment a plane circuit of sufficiently great extent carries an alternating 
current and the magnetic field due to currents thereby induced in the earth 
is mapped over the surrounding country. 

A method whereby the conductivity at any depth can be deduced from 
a knowledge of the induced field at the surface has been given by Slichter 
(1), but it is felt that this is not practicable and that it is likely that 
interpretation of the experimental results will depend on comparison with 
the fields calculated for some simple laws of variation of conductivity. 
Methods for computing the field are derived here for a uniformly conduct- 
ing earth when the circuit is so small that it can be regarded as a magnetic 
dipole. A circular current filament is then considered, attention being 
restricted to the vertical component of the induced field, this component 
being the one which would be measured in practice. 


2. The field induced by an oscillating magnetic dipole 

For a magnetic dipole whose moment M varies sinusoidally with period 
27/p, located at the point (0,0, h) outside a semi-infinite conductor of unit 
permeability and conductivity /, the inducing potential at the point whose 
polar coordinates are (p, 4, Z) is 
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where R? = p?+-(h—Z)?. 


{Quart. Journ. Mech. and Applied Math., Vol. IV, Pt. 1 (1951)] 
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t follows from equation (9.8) of (2) that since J,(Ap) satisfies equation 





6.11) there, the induced potential will be 


p_ ; 
Q2). M ; \ ‘ an Jy(Ap)A dx (¢ > 0), 
6+; 
here @? A? trikp c*. 
Writing /, w°=M, z W(h+), r vp, A xu, Where a? = 4zkp/c?, 
this becomes 
Q Mi > — af uJi(ur)ju du. (1) 
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Using the well-known result 


er DJy(Ap) da ({ <h), 


is can be expressed in the form 


M [ e-%*-OJ,(Ap)A dA. 


0 


When the oscillating dipole is at the surface we obtain for the magnetic 


eld components at the surface 


H rM,, ll ) = J,(ur)u2 du, (2) 
J V(r w)+u 

H VM, eal ) : J,(ur)u? du, (3) 
(ur r)+-u 


responding to an inducing field whose components at the surface are 


vl) 
H =) H, 0. 


“2 r 
? 


3. Evaluation of the integrals (2) and (3) 


Consider first the integral 


(w= b*) 


i 
if : J,(ur) du, 
(u- b?) +-U 
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where 6 is to be regarded as real and positive. We have successively Fo. 
f series 
l f¢ as as TT 
I 2 | {2u?-+-b?— 2u,(u?+-b?))J,(ur) du 
0 
l 9 tf The fi 
~ | Sf Ny2Lhe2 . 
pp | {,(u?+-b?)—u}udy(ur) du, comp 
0 induce 
: 1 referr 
since | J,(ur) du * indue 
0 | apped 
> ‘i The 
Using (u?+-b?)!—u b | e—"z J (bx) —, 
x 
0 
(3, p. 386), changing the order of integration, and integrating with respect 
to u, we find ‘ where 
2 c : : dx 
. 2 »2\—3 , 
I 1 | J, (ba) — f(a2+-72)-9 
i b * Cx x 
0 
l 21d fF J,(be 
et | -— dx 
r brdr J (x?+r?) ; 
: using 
1 2ld : | 
+ (1, (br) K,(4br)} 
r brdr*’* , 
(3, p. 435), which can be written 
1 21d /fl—e™ 
I —}. (4) | 
r brdr br 
The restriction that 5 is real can now be removed by analytic continuation, 
provided only that its real part is positive. 
Finally, writing Bessel’s equation in the form —_ 
¥ ldj{d | 
u>J,(ur) r—J,(ur)|, 
rdr\ dr | 
we can derive H, from the integral J, giving wher 
ldf{ d{l, 2 d[1—e-\) a1 = 
H. M, r — ; (9) 1] 
. rdr| dr\r- brdr r an 


where b = 7, that root being taken for which the real part is positive. 
Performing the differentiations we can also write 
1 18 . 2be "( 4. 9 9 \) the ft 
Lo 


; (6) 
| , are té 


or | (br)? T (br 


H = —clf, 


ee nn 
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ly For small values of r (up to unity, say) it is more convenient to use the 


series eX pansion 





H 21M, b* ~~ yn +1)*(n 3) (py 7 (7) 


a (n-+-4)! 
The form of H, given by (6) shows that when r is large enough the vertical 
mponent of the induced field is effectively equal and opposite to the 
nducing field. This is likely to be of importance in the experiment 
eferred to in section 1 since it is true irrespective of the shape of the 
nducing circuit (which can be regarded as a dipole distribution) and also 
ppears to be true for the case of variable conductivity. 
The evaluation of H, proceeds in a similar manner. We have 





cM 
H rM,—, 
or 
espect 
(w?-+-b?) l 
vhere i) : : ud,(ur) du 
J (a b*)*+-u 
‘ dx f 
9 e~" J,(bx) ud, (ur) du, 
x 
0 0 
ing (3), p. 386. Therefore, 
i l > > dx 
ri) 2 J,(bx) f(a2+-r2)-2} 
F ‘“ dx r 
2d £ Jit 
{ 2\ wa dx 
r dr J (x*+7r?)! 
1ation, | 2d : 
— {1,(46r) kK, (46r)} 
f ar - a 
sing (3), p. 435. Hence 
dfld ; 
H. rM,, 2 — {L,(4b6r) Ky (br)} (3) 
dr\|rdr**** i 
nere h z=. 
Performing the differentiations and using the recurrence relations for J 
nd K (3, p. 79) we can also write 
ive 6 


(Lh Ko—L, Ky) ’ 


the functions involved being Kelvin’s ‘ber’ and ‘bei’ functions, etc., which 
re tabulated. 


b 
- (9) 


l 1h : : 
H My IF K, 9 (Lo Ay I, Kk) —— 
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For very small values of r it is more convenient to use a series expansion | It fo! 
in ascending powers of r, which can be derived from the result given by 
Watson (3, p. 441), 


hor 


9 Pf 
I,(z)K,(z) = — | K,(2zsin @)cos 26 dé, 


by expanding K,(2zsin @) and integrating term by term. The same series 
is also obtained by putting » = 0 in (12) (see later). 

The values in Table I were computed by means of the above formulae, 
the factor aM, being omitted throughout. 


TABLE I 





using 
Induced field 
Inducing —— 
field 
r H, a. H, 
o*2 125°0 0°1646+1°06252 | 0°0313+1°24032 W hei 
°° 15°625 0°1429+ 0°44021 0°0455+ 0°60611 
06 4°630 0°1231-+0°23611 0°0534+ 0°38941 
08 1°953 0°1054+0°13721 0°0576+0°27771 
I°o 1000 0°0896 + 0:08082 00593 + 0°20571 
1°5 0*296 0°0576+-0°01491 00569 + 0°11 301 
20 O°125 0°0347—0°008 32 0°0497—0°06417 
4. Evaluation of the integral (1) 
Consider the integral 
x 
% du 
I e-"= J, (ur) . (b real). 
1 0 ) 9 : 
. (u? +6? and 
0 
This can be written 
7 x 
ey r du 
9 —U(z+4 ir cos d) 
I, | dd | , —_ 
TY. J \(u*+-6*) 
0 0 
7 x x 
a i ‘ Se 
=- | dd | e—uc+ire Os d) du e-un, (Dx) dx. 
7. . 7 whe 
0 0 0 


west ( dd f tbe, 


wt-+K 


Ix, 


e © 


0 0 
where x = z+ ircos¢. 
The second integral is known (3, p. 436, (7)), so that 
I, = 3 [ {Hy(bx)—Yo(bx)} dd. 


0 
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b2 
iF 
here 
] 
. z a 
F(x) - a — 
Lica i Ue 1 3° 
l l 
; Jj 
Dre i 
] l 
( 
: t 3 
( 0 5 log 
ind y is Euler’s constant. Now 
a" P. (pL) 
R (n DP (yu) 
where R2 2+-7?2 and pu z/r, 
R. > _ R3. 2 
I, p34 P, psi F. 
; 3 37.5 
) 
P, ~*~ (Cc. 
R2°24, 
< ( ] "2 2 
hay DOM lin'(m 1)! 
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b*)te-“= J,(ur) du (10) 


Y,(bx)} dd (11) 


This can be written 


{b(z+-ir cos d)| dd, 
=o, T 
l m ; (4a)2™ 

x)logz+— > (—1)"C,, —2 , 
a 7% “m!(m-+1)! 
| . —y-+log 2, 

m ~ 2(m-+1) 
» 

| (z+ir cos ¢)”" dd, 

0 
U fee 

(z+ir cos d)-™+» dd, 


so that the integral (10) becomes 


, &.F, 
SS goed 
f2m+2 2m A 
log b)( ])™ a R: I 2m _ 
22™m!(m- 1)! 
lf 
(z+-1r cos ¢)?” log(z+-tr cos d) dd. 








A. N. GORDON 





112 


Since both sides still converge when b is complex, we can remove the 
restriction that 6 is real. 

The function 
a > i , ‘ d 
i (z+7ircos¢)" log(z+-ir cos ¢) dd = 


7 dn 


(RP, (H)5, 


is a special solution of Laplace’s equation (4, p. 172). It has the value 


+z 2n—1 2n—3 > 
R*P (u)log dt a cece Pot (ie | 


2 \ 1.2n 2(2n—1) n(n-+-1)} 


The required potential (1), omitting the factor M), is thus 


Q, . {2u?+-b?—2u(u?+-b?)!e-"=J(ur)u du 
52 
0 
12P, 2a 
Plt) 5 a = (u?+-b?)'e-“" J, (ur) du. 
>2 2 2 B22 
0 
Now ‘ 
(RP, (u)} = rR" P, 1(u), 
Cz 
-,{R (n DP (u)} (n | 1) R-+2)P) os (p); 
and —{F,} nF,_,+ RP, _1(u). 
Hence, finally 
is 4 alias ai 
Q, = 2 BRP, ORF, +. if og 
13.5 3°.5.7 Te Ce ce - 
ae 9 9 
S (C,—logb)(—1)™ a 1) b2m RIm-2P 
—, 22"™m!(m-+1)! a 
m 


Se (= 1mm ' ~~ 
bi >. 22mm! (m + ])! (2m(2m 1) Fym—2+ (4m- 1) vem von 2s: (12) 
1 


m 
This formula is suitable for the calculation of the fields for values of R not 
substantially greater than 2. 

Proceeding in the same manner and using the asymptotic expansion for 
H,(z)—Y,(z) (3, p. 333 (2)), we arrive at the formula 


i~ 


2 P, _ 12P; , 9 > (—1)" (2n)!(2n-+-2)! Pons 


OR BRE” ZL, (2n—1)2?"n! nm! BARI” 
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which is suitable for values of R not substantially less than 10. In order | whe 
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» the | to compute the field components for intermediate values of R numerical 
integration must be performed from u = 0 to u (= 2, say), the integral 
from x to 00 being expanded into a series of integrals of the form 

r du , ¥ du 
| Jy(ur)e-* —, K, = | J,(ur)e-“ 


3 u* j u® 


I zx 


These can be linked with J, and J,, by means of reduction formulae, where 


ue 
J, and J, are given by 
| iy l a yn +1 
I, = | Jo(urje-“ du aw te RP, (p), 
R — (n- 1)! 
be n=0 
i du 
I, J, (ure - 
: u 
a gn 
log 2—y—log x—log(z+R)— Y (—1)"—, R*P,(u). 
wy n.N: 
n=1 
5. The circular current filament 
The simple form of H, in equation (5) enables the same component to 
be calculated in the case of a circular wire of radius d which carries an 
alternating current cJ. For such a circuit, whose centre is at (0,0,h) and 
whose plane is parallel to the bounding face of the conductor, the inducing 
| potential is 
Q) 2aTd | eAh-9J (Ad)Jj(Ap) dA (C <h). 
0 
In terms of the reduced units introduced in section 2, and with a = ad, 
the induced potential is found to be 
(a*-+ b? u = 
Q) 2rla : ) e~“=J,(ua)J,(ur) du, 
, = b?) + U 
J Ve 
0 
where b? i, so that 
eQ), * /(a?+-b?)—u : 
(12 H 27 laa ——- e-“=J,(ua)J,(ur)u du; 
. o¢ :  (a* ; b*) +- U 
0 
Rnot this can also be expressed in the form 
é f ,(a*+b?)—u . 
ion 101 H 27laa 5 5 e~2J,(ua)J,(ur) du 
. ca J .(a*+-b*)+u 
U 
on : 9 9 
‘ (a*- b? u ” 
Taxa dd \ > . e~"* J,(uR) du, 
ca . (a? +-6?)+-u 
0 0 
1 orde! where R2 a*—2ar cos¢+r? 
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In the special case z = 0 we have 


on 


H. te | dd iF A " Jo(uR) du 
a oa (w?-+-b?)+-u 
0 
af (1.21 d /l—eR\) 
tia. | we Id, I 
ea | \R* 6 RaR| bR )}' ie 


0 
as in (4), 
Expanding the integrand in ascending powers of R and performing the 
differentiation with respect to a, 


_ n(n-+-2) Lees oe 
oe Ma ” (n--3)! be (@— 17) X,-2 +X}, (14) 
where X,, | R” dd. 
0 


These integrals can be computed from the tabulated elliptic integrals 
X_, and X,, using the results 


X_ = Qn, X, = 2n(r?+a?), 
together with the reduction formula 


(n+2)X.5 = 2(n+1)(r?+a?)X,, —n(r?—a?)2X,,_». 
The expression (14) is suitable for computations involving small values 
of r and a. For larger values an alternative expansion can be derived as 


follows. We write (13) in the equivalent form 


9 
27 





l 2 2 d{ /2 K,(bR)\ 
H. Ia i iets d 
, oa | R bR> dR Rs 7 DR)! | ? 
rfl 2 /2 K,(bR) 
a Oe 7 1 QF d 
a) [RO BR Te (bP) | $. 
To evaluate this integral we note that 
K,(bR) 2\s—.,, 3\(7 
(bR)! \ (jae 3 (2n+-3)C,,(cos 4)K,, ,;(ba)L,,.,(br) (a > 7), 
where (1—2h cosd+h?)-? ¥ hC,,(cos 4), 
n=0 


(3, p. 365, 


8), and interchange of J and K should be made when a <r. 
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The value of , 


I, ( C,, (cos ¢) dd 


0 


which is required, can be obtained from the recurrence relation 
27 
[,—I,-2 = (2n+1) | P,(cos 4) dd, 
0 


{1.3...(2n—1)}* 
92ny 1. ) (2n-+-1)2z7, n even, 


22"! n! 


I 0, n odd. 


x the 


(14) 
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ELECTROMAGNETIC INDUCTION IN A UNIFORM 
SEMI-INFINITE CONDUCTOR 
By A. N. GORDON (Imperial College, London) 


[Received 30 May 1950] 


SUMMARY 
A systematic treatment of the problem of the induction of currents in a uniform 
semi-infinite conductor by external magnetic or electric fields is given, attention 
being drawn to the complementary nature of the magnetic and electric cases. 
A one-dimensional heat flow analogue is then derived which enables the familiar 
methods developed in connexion with heat conduction to be applied directly to 
the solution of corresponding problems in electromagnetic induction. 


1. Introduction 
Tuts paper should be regarded as a sequel to a recent contribution by 
A. T. Price (1). Since the method of approach is rather different, however, 
it has been considered advisable to develop the subject afresh, without 
reference to Price’s work. 

With the usual notation, when displacement currents are neglected, 
Maxwell’s equations for the electromagnetic field inside a uniform source- 
free medium of conductivity « and permeability pu are 


curl E - _ (1.1) 

c at 
curl B - — E, (1.2) 

. 

div E — 0, (1.3) 
div B = 0. (1.4) 

Elimination of E yields for B the equation 
-SS=. (1.5) 

c* a 


while elimination of B leads to the same equation for E. 


2. Solution of (1.5) subject to (1.4) 


In rectangular coordinates, the most general solution of (1.5) is 
B = iW,+-jW.+kwW,, (2.1) 
where i, j, kK are unit vectors along the axes and W,, Wy, W, are arbitrary 
solutions of the scalar equation 
ae Low 
VW = -—., 
a ot 


(Quart. Journ. Mech. and Applied Math., Vol. IV, Pt. 1 (1951)] 
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ELECTROMAGNETIC INDUCTION 


where x = c?/4rpK. (2.3) 
The additional requirement (1.4) imposes the condition 
C W, om, Mn 0). (2.4) 


Ox oy Oz 


so that only two of the scalars W,, W,, W, are independent. Since there is 
no loss of generality in replacing W, W, by eW,/éz, eW,/éz, it follows from 


(2.4) that 
W. if ow, , om 
3 dx oy ° 


Strictly speaking an arbitrary function of x, y, and ¢t should be included 
in W,, but this special type of solution, which leads to fields independent 
of z, will not be considered here. Equation (2.1) then becomes 


™ (ie Kt) - (1% th 
Oz Ox Oz oy 


curl(jW,)+-curl(iW,). 

Since W, and W, are independent, it follows that the most general 
solution can be obtained by combining simple solutions of the form 
eurl(iW) and curl(jW), to which may be added, by symmetry, curl(kW). 
However, using the identity 


curl grad W curl iT rv : +e} - 0, 


| 


Ox Cy Oz 


any one of these types can be expressed in terms of the other two. Further, 


curl curl(kW) curl 4°" +12), 


Ox oy 


since 


two independent solutions are 


eurl(kW); curleurl(kW), (: 


bo 
or 
— 


the second of which can also be written 


grad|‘ =) : We 
Oz x Ot 
in virtue of (1.5). 

Equations (1.1) and (1.2) indicate that if E is identified with either of 
the solutions (2.5) then B corresponds to the other. This reciprocity is 
merely the degenerate case of the same feature in the solution of the 
electromagnetic wave equation (see, for example, Stratton (2), pp. 27, 28). 


3. Application to the semi-infinite conductor 
Suppose that the region z > 0 is occupied by conducting material while 


an external magnetic inducing system is confined to the region z < —A, 
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where h is positive. The equations to be satisfied by the magnetic field are 


V?B = . = (z > 0), 

x ot 
V5B = 0 (—-h <z< 0), 
div B = 0 (—h <z < 00), 


subject to the usual boundary conditions that B., H,, and H, are con- 
tinuous across z 0. 

When a suitable solution has been obtained, the electric field inside the 
conductor follows uniquely from (1.2). Outside the conductor we must use 
(1.1), which leaves E indeterminate to the extent of the gradient of some 
scalar potential function. 


To solve the equations we must first find a suitable solution of (2.2). If 


we seek solutions of the form 
W Z(z,t)P(x,y) (z > 9), (3.1) 
Z,(z,t)P(x,y) (z < 9), 


the boundary conditions (which apply over z 0 for all values of ¢) will 
involve only the functions Z and Z,, and substitution of (3.1) into (2.2) 
indicates that 


eP #@P .. 
af. ———_. -. {8 P 0. (3.2) 
ox" cy” 
CZ is 1oZ 
— — A?Z , (3.3) 
G2 x ct 


where A is a positive parameter. 
Outside the conductor corresponding solutions of the same form are 
obtained by putting x = 0 (i.e. a = 0), giving 


W = CetP(x, y), (3.4) 


where C is a function of ¢ only. 


4. Solutions of magnetic type 
Taking B = curlcurl(kW), where W is defined by (3.1) or (3.4), we get 
B (‘ ZOoP @Zal 


102 O2 C2 cy 


C ) > . y > g > \ 
B= Ae | ie ‘3 AP, Bee AP 


ox CY 


22?) (z > 0), (4.1) 


grad(Ae-“4 Be*)P (z < 0), (4.2) 


where A and B are functions of ¢ only. 
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Corresponding values of the electric field, derived from (1.2) and (1.1) 


respect iv ely , are 


E fin: wail 50) (z > 0), (4.3) 
c ot \cy Ox j 
1 {dA dB oP oP 
> e> ,O)4grad (2 <0), (4.4 
cA\ dt it Ney Cx ready ( \ ) 


where % is an arbitrary potential function. 
The boundary conditions for 1/,, H,, B, at z 0 yield 


{(t)+- B(t) =) . (4.5) 
pb Cz} o 
1(¢)— Bit) \(Z)o (4.6) 
and. eliminating A(t). 
1/cZ . ‘i 
—| A(Z)y = 2A(t). (4.7) 


Examination of (4.2) indicates that the term with exponent (—Az), 
which tends to zero as z tends to infinity, may be interpreted as an 
nducing magnetic field, the term with exponent Az being the induced 
field. For a given problem, where A(t) is a known function of t, (4.7), 
together with the physical requirement that Z tends to zero as z tends to 
infinity, may be regarded as leading to a unique solution of (3.3). It also 
ollows that (4.6) ensures the continuity of E,, E, at z = 0, provided that 
there grads has at most a normal component and hence corresponds to 
some electric charge distribution in equilibrium on the boundary. Apart 
from this case, which corresponds to a superimposed electrostatic problem, 

can be taken as zero. 

In order to solve a given induction problem by the method described 
ibove it is first necessary to express the inducing potential in the form of 
1 sum (or integral) of expressions of the type 

Ce-“ P(x, y). 
This can always be done, in principle, at least, by means of Fourier analysis. 
For example, the potential of a unit pole can be expressed in this form by 
means of the well-known result 
ee : | e-=J,(Ar) da. 
0 
5. Solutions of electric type 

In exactly the same way, starting with E = curlcurl(kW), 

E obi 


Cz Ca 2 oY 


, BLP) (z > 0), (5.1) 


E= A ie AP) Be SAP) 


Cx CY ox Cy 


orad{ (Ae Be*)P} (2 <0), (5.2) 
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the corresponding values for the magnetic field being 


B= “ae ae, 0) (z > 0), (5.3) 
x \oy Ox 
l dA Hans dB y-\ (eP oP e 
= — 71 —_@W »,—-+~,O0}]+grady (z < ( dD. 
” al di* ' dt. le Ox ) er 7 


where ¢ is an arbitrary potential function. 
In deriving (5.4), in order to maintain the symmetry between the two 
types of solution, the accurate equation 
1 cE 


c ot 


curl B 


has been used in place of curl B = 0. 


The boundary conditions for both E and B are all satisfied if we take — 


% = 0 and 


A(t)-+ B(t) (<) (5.5) 
CZ So 

dA dB - ie 

aS 4AK(Z)p. (5.6) 


This second type of solution is seen to correspond to the induction of 
currents by an external electric field. Since « is large, equation (5.6) 
indicates that (Z), is very small and (5.5) then indicates that under 


ordinary circumstances A(t) ~ — B(t). 
The associated electric surface charge is given by 
do , " ; aig 
7 -1, = —KE, -Kr?(Z)y P, (5.7) 
C 


where 7, and FE, are the normal components of the current and field inside 
the conductor at the surface z = 0. 


6. Reduction to a problem in heat conduction 

The formal resemblance of equation (1.5) to that governing heat con- 
duction has often been noted and, in some simple cases, known solutions 
of heat conduction problems have been applied to those in electromagnetic 
induction, but in general the boundary conditions of an induction problem 
have no counterpart in heat conduction. For the type of solution con- 
sidered in section 4, however, it will now be shown that a complete 
analogy does exist. 

According to section 4, the currents in a conductor due to an inducing 
field which is derivable from a magnetic potential 


Q = —Ae-“P(a, y), (6.1) 
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can be expressed in terms of a function Z(z,t) which satisfies 


SY Mme Se, (6.2) 
0z* x ot 
while at the boundary z 0. 
eZ - - 
pAZ 2A (t). (6.3) 


Cn 


There is the further requirement that Z should tend to zero as z tends to 
infinity. 


Making the substitution Z Ye-"at, (6.4) 
a oy 1 oY " 
we have to satisfy ; - 2 : (6.5) 
oz x Ot 
. oY , 2 af _ a 
with pad 2pue*™ A (t), (6.6) 
C2 


at 2 Q. 

Equations (6.5) and (6.6) have the same form as the corresponding ones 
in problems of heat flow in one dimension in a semi-infinite conductor, where 
the temperature 7’ satisfies the equation 

eT 1¢éT 
oz*@ ok Ot’ 
« being the diffusivity. 

For radiation at the surface into a medium at temperature f(t), there 

is the surface condition 
ol 


Oz 


-hT = —hf(t), 


where / is a constant related to the emissivity of the conductor. Solutions 
of a variety of problems of this type are known (see, for example, Carslaw 
and Jaeger (3)). It may also be noted, in passing, that the general question 
of heat conduction in a semi-infinite conductor (i.e. not limited to one- 
dimensional flow) can be treated by the methods described above. 
A slightly different procedure will now be adopted which has the 
advantage that it can be modified to cover the electric case as well. 
Vonsider the function U’ defined by 
a (6.7) 

Oz 
This satisfies the same equation as Y, namely (6.5), throughout the 
conductor and its value is given over the boundary by (6.6). We can 
therefore identify U with the temperature of a conductor whose surface 


temperature is prescribed. 
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If we assume that there are no residual currents in the conductor other 
than those due to (6.1) we can now write down immediately 


— 


Ce N77 
0 


oY — 2 — ee , 
— pA P Al 4 | VAMt—(27/4axu*)} p—u® ay (6.8) 
J xu~ 


which is seen to satisfy all the requirements of the problem (cf. (3), p. 44, 
equation (1). The lower limit there is modified, however, since A(t) is not 
necessarily zero prior to t = 0.) 

An example of some interest is the sudden creation of a magnetic field 
outside the conductor. Taking Heaviside’s unit function H(t) for A(t), 
(6.8) becomes, when use is made of (6.4), 

= — paAZ 4H [ ¢-w0rs iw’) du, where B = 32(at)- 


C2 N77 
. 
R 
p 


(6.9) 


The evaluation of the current given by (4.3) is quite simple since it 
involves @Z/ét only, and differentiation of (6.9) with respect to ¢ gives for 
eZ ct the linear first-order equation 


- » ‘eis \3 
02\ Oz C KC 


\ 


7. Free decay of a ‘magnetic’ current system 

In order to deal with the decay of a magnetic system of currents under 
the action of its own field it is necessary to put A(t) = 0 in (6.6). We 
require then the solution of 


ey 1¢6Y 
= ; (7.1) 
C2" x ot 
, ; oY . m 
for which — — pA} 0, (7.2) 
oz 
at z = 0. The heat analogue will be the decay of temperature in a semi- 


infinite conductor whose plane boundary is maintained at zero tempera- 
ture. 

The general solution of (7.1) subject to (7.2) is well known and is given 
by 

. eY 


C2 


- prY | F(d)sin (fz)e-#°™ dd, 
é 


where F(¢) is arbitrary. Using (6.4) this becomes 


on — pA | F(d)sin (dz)e“#? + dd. (7.3) 
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ther Solving for Z, we obtain 
. .. b> 4+-pAsin dz 
? 9 © COS O2 LA SIN Oz 2,2 - 
Z k (hd : ms ¢ (2 +A*)at dd, (7.4) 
ia ge My? 
0 


the complementary function (which involves e4#*) being rejected since it 





. 44 ‘ : , _ . . 
nereases indefinitely with increasing z. The free modes, defined by 
Z (fd cos z+ pA sin dz) (p?+A*)att | (7.5) 
field 
114 re immediately evident. 
: The free decay of a current system, which is everywhere defined at 
0. can now be investigated. Let the value of Z(z, t) at ¢ 0 be denoted 
y Z,. Then, according to (7.3), F(¢) is to be determined from the integral 
6.9 . 
juation 
dZ , ; : 
; S— nrAZ, | F(¢)sin dz dd. 
ce it dz 


‘S fol 0 
Inverting this by Fourier’s integral theorem we get 


6.10 


F(d) = — | — pAzZ)sin dz 
nder By . 
We : (A cos dz+-pAsin ¢z)Z, dz, 
(7) rovided Z vanishes at z 0. 
If, instead of Z, €Z et is given at t = 0, we can proceed in exactly the 
te same way. 
(7.2 
It is instructive here to consider again the case of the sudden creation 
emi- | ofa magnetic field outside the conductor. It is tacitly assumed in neglect- 
pera ng displacement current that the time rate of change of the fields is 
| sufficiently small, since it is the term involving éE/ét that is neglected in 
riven juation (1.2). Sudden creation of fields violates this initial assumption 
nd the sudden appearance of surface currents in the magnetic case and 
surface charges in the electric case is the price we have to pay. Once these 
| surface distributions have come into being, however, they ensure that 
subsequent changes proceed sufficiently slowly for the initial assumptions 
to be valid. It should also be noted in this connexion that the parameter ¢ 
7.5) must not be too large and that in integrals like (7.4) we rely on 
(7.3) their rapid convergence to make contributions from large values of ¢ 


elatively unimportant. 
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Consider then the inducing field 


H = an rel -P) (t > 0), 

ox” Oy 
suddenly applied at ¢ = 0. It is well known that under these circumstances 
a surface current is instantaneously set up on the face of the conductor so 
as to completely shield the inside of the conductor from the applied field. 
The density per unit area of this current is given by 


oe? a ees 


x F y y : x 
2a . 27 


A much clearer picture of this process is obtained if the external field 
is regarded as approaching the conductor with a very large velocity in the 
form of a travelling discontinuity, on either side of which the magnetic 
field can be represented as the gradient of some scalar. Equation (1.1) then 
indicates that the electric field is everywhere zero except at the discon- 
tinuity, where it is infinite. It is this infinite electric field which produces 
the surface current on the conductor. Thereafter we can regard the field 
as completely reflected so as to set up the initial induced field 

_ oe ail a?). 
Ox oy 


The surface current now diffuses into the conductor and we have a case 








| 


; 
' 
| 


of free decay. When the current has everywhere decayed to zero, however, | 


there remains a stationary induced magnetic field 
H = ] —F ode oP ep XP 
1-+p Ox Oy 
outside the conductor, and a field 
en wl 
l+-p Oa’ dy’ 


inside, so that ultimately 


Z — _4 {yy ‘. 
Ap+l 7 

Z ‘ 
and d —ApZ = Que-™, 

dz 
> ‘ oZ ‘a a - ; (62 at or 
Putting ——AuZ = Que-*+ | F(¢)sin (dz)e (f2-+A2) dd, (7.6) 

oz 


Pa 


and noting that at ¢ = 0, Z vanishes identically throughout the conductor, 


we have 


i F(¢)sin dz db = —2pe-™, 


0 
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of which the solution is 


. 4 
I (d) P E _ 
a ¢?+)? 
Substituting this into (7.6) and solving for Z gives finally 
Z 1 2 ore 4H [ dcos(gz)+Ausin($2) » igrsrnaig de, 
Ap+l 3 (f2+-A?)(d?4 Au”) 


0 
If (7.6) is differentiated with respect to t, the equation (6.10) is again 
btained after some reduction. 
It has now been shown in this section that the decay of a given current 
system can be described in terms of a function Z(z,t) whose initial value 


Z, is everywhere given. If 


se AuZ, Q, 
the process is straightforward since the magnetic field as well as the 
current system decays completely. Such a state of affairs could be brought 
about either by internal causes or by an external magnetic field which 
ceases to change as it passes through the value zero. 


If. on the other hand, dZ, 


L AuZ, ~ 0, 
there is a permanent residual magnetic field which must be taken into 
vecount in calculating the decay of the current system. This state of affairs 
can only correspond to an external field which ceases to change at ¢t = 0, 
maintaining its non-zero value subsequent to t= 0. Sudden changes in 
the external field and the accompanying surface charge can be treated 
either by the methods of this section or as special cases in section 6. 

Finally it will be pointed out that the free modes given by (7.5) also 
apply to a conducting slab occupying the region 0 < z < 2h. In place of 
an unrestricted continuous range of values, however, ¢ is limited to the 
discrete set defined by étan(¢h) at 
8. Induction of currents by external electric fields 

Elimination of B(t) from equations (5.5) and (5.6) indicates that, for the 
electric case, solutions of (3.3) are required subject to 


“Z r, dA 
4nkAZ 2 (8.1) 
ozot dt 
at the boundary z = 0. The surface charge given by (5.7) becomes, with 


the aid of (8.1), 


2 | 


AP 1/eZ 
4 (Z) . (8.2) 
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The method of section 6 leads to 


x 


n97 9 
eZ . 4u , = Si a 
—— —4rkAZ f [ 4 (/ — = e-APz*/dau*)—* Jay 
oz0et Vi. 40u* 
0 

where the dash denotes differentiation with respect to 7. 

Since « is large even for only moderately good conductors, for external 
fields which do not vary too rapidly, to a high degree of accuracy this 


reduces to 


Z _ [ a'( m le (APz*/Aocee®)—e® Gy (8.3) 


4au? 


with o= —Ai(i). (8.4) 


Thus the electric field inside the conductor is very small, but this is not 
true of the current. 

Equation (8.3) is unsuitable for dealing with sudden changes in the 
external electric field, which require special treatment. The remarks of 
section 7 in connexion with sudden changes in external magnetic fields are 
also applicable here. Introducing again the idea of a travelling electric 
field, the associated magnetic field is everywhere zero except at the electric 
wavefront, in which it lies. The penetration of the electric field alone into 
the conductor can only set up a finite current system and cannot lead to 
an instantaneous surface charge. On the other hand, in the infinite mag- 
netic field carried along by the electric wavefront, we have an agent 
capable of producing the infinite current ‘loops’ at the surface necessary 
to bring such a charge into existence, and the arrival of the electric field 
is accompanied instantaneously by just that surface charge which prevents 
further penetration of the fields into the conductor. It is important to 
notice that the redistribution of surface charge which takes place is not 
due to a surface current, but to a system of current ‘loops’ confined to the 
surface layer of the conductor, the flow being normal at the surface. As 
such, this behaviour differs from that shown by equation (8.3), where there 
is appreciable current flow at a finite depth below the boundary. 

Up to this point the electric and magnetic cases are completely comple- 
mentary. In the magnetic case, however, the initial surface current diffuses 
into the interior of the conductor and, in due course, the incident magnetic 
field is able to penetrate. In the electric case, on the other hand, the 
surface charge remains, since a conductor cannot sustain an internal 
stationary electric field. 
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9, Free decay of an ‘electric’ current system 


Putting dA/dt = 0 in (8.1), we find that the equation governing free 


decay 1s e2y l oy 
°° ans ? (9.1) 
0z* x ot 
eZ 2 
subject to 4nkAZ = 0, (9.2) 
ternal — 
v this at the boundary z = 0. 
Proceeding as in section 7, the general solution can be written 
9 eZ ‘ eo , — 
(9.9) 4mkAZ F(¢)sin(dz)e“P* +"! dd, (9.3) 
Cx 
Uv 
(8 4) , * 4A sin(dz) — Bd cos(dz _— 
8.4 or Z, F(¢) me )—f ; * (fp? +A*)at dd, 
: B®? +-(47rKA)? 
0 
IS not . a 
where B = (¢?-+-A?)a. 
The free modes of decay defined by 
n th 
lL f y ° ‘db 2 
rks of Z sin(dz) — B cos(d)} pr 
ds are 4akr 
lectri¢ are immediately evident. 
lectric To the same order of accuracy as was adopted in the magnetic case, the 
e into | situation is much simplified, since a solution of (9.1) is required which 
ad to vanishes at z 0 (i.e. only the term 47xAZ is retained in (9.2)). This can 
mag either be written down as a simple modification of (9.3), or can be exhibited 
agent in the equivalent form ((3), pp. 34 and 40), 
‘ssa;r\ 
» fie , l — . 
e field } F(z+-2( vt) Pue u? du. 
‘vents Var 
int to s . : . A — , 
where F(z) is an odd function of z. It is easy to verify that Y reduces to 
is not : : - 
F(z) when ¢ 0. 
to tne A Z . . ° 7 ° 
It follows that any choice of an odd function F(z) will lead to some 
\ * 
Cc. | » ‘ . 
freely decaying system; for example, 
ther . ; ; 
F(z) = sin(¢z) 
mple gives Z e(P? +A*at sin(dhz), 
ffuses ; , 
which corresponds to the free modes of decay. 
neti . 
1, the Another example is F(z) ze—mpxdst|c* 
ternal which leads to 
Z (1 bt)—!eVetsAmpKye mK pLb2*/(1+bl) 
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At any time ¢, Z has a maximum at a depth given by 
_ _ {c%(1+-bt))4 
\ | 2apuxb | 
which maximum moves progressively into the conductor. 

The magnetic and electric cases have so far been considered separately, 
By expressing an arbitrary non-divergent current system in the form of 
a triple Fourier integral, it is not difficult to decompose it into constituent 
electric and magnetic parts, so that, in principle at least, its decay can be 


described by the above methods. Ifthe separate parts do not automatically 


satisfy the necessary relations at the surface, the work will of necessity be 
complicated by possible surface charges and currents and permanent 
external fields. 


10. Conclusion 

Inducing fields that have already received fairly detailed treatment are 
those produced by (a) an infinite straight alternating current filament 
(Price (1)) and (5) an alternating magnetic dipole and a circular alternating 
current filament (Gordon (4)). In these, however, attention is mainly 
confined to the external induced field (the latter can be determined by 
means of (4.6) or (5.6) once Z(z,t) has been found). 

Another class of problems that can be solved fairly easily involves 
isolated magnetic poles or electric charges in uniform motion perpendicular 
or parallel to the surface of the conductor. 
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